NPS  ARCHIVE 
1963 

DEFIBAUGH,  C. 


LIBRARY 

US.  naval  postgraduate  school 

MONTEREY,  CALIFORNIA 


INELASTIC  COLUMNS  SUBJECTED  TO 


PULSE  LOADING 


Carl  P.  Defibaugh,  Jr 


INELASTIC  COLUMNS  SUBJECTED  TO 


PULSE  LOADING 


by 

Carl  F.  Defibaugh,  Jr0 
Lieutenant,  United  States  Navy 


Submitted  in  partial  fulfillment  of 
the  requirements  for  the  degree  of 

MASTER  OF  SCIENCE 
IN 

MECHANICAL  ENGINEERING 

United  States  Naval  Postgraduate  School 
Monterey,  California 


19  6  3 


Library 

US.  NAVAL  PO'V.  RAOUAIE  SCHOOL 
M.-NiiU.LY  CA'.i  -ORNIA 


INELASTIC  COLUMNS  SUBJECTED  TO 
PULSE  LOADING 
by 

Carl  F.  Deflbaugh,  Jr. 


This  work  Is  accepted  as  fulfilling 
the  thesis  requirements  for  the  degree  of 
MASTER  OF  SCIENCE 
IN 

MECHANICAL  ENGINEERING 
from  the 

United  States  Naval  Postgraduate  School 


ABSTRACT 


This  paper  reports  an  investigation  of  the  dynamic  behavior 
of  a  pin-ended  Inelastic  column  subjected  to  a  half-sine  dis¬ 
placement  pulse  at  one  end  applied  with  small  eccentricity. 

The  rates  of  end  displacement  are  such  that  the  effects  of 
axial  inertia  are  considered  negligible.  The  column  is  re¬ 
placed  by  a  lumped  parameter  mathematical  model  and  the  equa¬ 
tions  for  the  model  are  solved  by  a  high  speed  digital  com¬ 
puter.  The  stress-strain  relation  during  loading  is  assumed 
to  be  of  the  Ramberg-Osgood  type.  Unloading  occurs  along  a 
line  with  slope  equal  to  that  of  the  tangent  at  the  origin. 
Slenderness  ratios  of  50,  70,  and  100  and  two  different  values 
of  eccentricity  are  considered  in  the  investigation.  The  max¬ 
imum  load  supported  by  the  column  and  the  residual  lateral 
deflection  at  the  center  of  the  column  are  found  as  functions 
of  the  displacement  pulse  height  and  of  the  displacement  pulse 
duration. 

It  is  shown  that  an  inelastic  column,  like  the  elastic 
column,  can  support  dynamic  loads  considerably  in  excess  of 
the  static  load  capacity.  It  is  also  shown  that  within  the 
range  of  pulse  durations  studied,  the  smallest  residual  de¬ 
flections  generally  accompany  the  largest  maximum  loads. 
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INTRODUCTION 


In  the  study  of  structures  the  column  ranks  as  one  of 
the  more  provoking;  elements  as  judged  from  the  amount  of 
literature  published  on  its  behavior.  The  development  of 
elastic  column  theory  is  comparatively  advanced  and  dates 
back  to  1757  when  2uler  published  his  analysis  of  the  strength 
of  columns.  In  more  recent  years  Inelastic  column  theory  has 
been  developed  to  the  point  where  a  satisfactory  explanation 
has  been  obtained  for  the  phenomenon.  Shanley  /l/*  has  been 
credited  with  revising  and  reinterpreting  inelastic  column 
theory  from  previously  accepted  erroneous  conjectures.  Util¬ 
izing  studies  of  a  highly  idealized  column  he  concluded  that 
th9  tangent  modulus  load  is  the  maximum  load  at  which  an  in¬ 
itially  straight,  centrally  loaded  column  will  remain  straight; 
loading  beyond  this  load  causes  bending  which  produces  per¬ 
manent  deformation. 

Tne  dynamic  behavior  of  structures  has  become  an  increas¬ 
ingly  important  engineering  consideration  in  recent  years  and, 
with  the  aid  of  more  sophisticated  devices  for  the  problem 
solution,  many  more  phenomena  can  be  Included  in  an  analysis 
than  could  previously  be  reasonably  considered.  The  behavior 
of  dynamically  loaded  elastic  columns  has  received  consider¬ 
able  attention  and  a  number  of  investigations  have  resulted, 

Hoff  and  associates /2, 3/  considered  the  case  of  an  elas¬ 
tic  column  of  uniform  cross  section,  initially  curved  in  the 

*  Numbers  in  /  /  refer  to  entries  in  the  bioliography 
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shape  of  a  half-sine  wave,  subjected  to  constant  velocity 
loading  such  as  that  encountered  in  compression  tests  in 
commercial  testing  machines.  It  was  shown  that  the  maximum 
loads  supported  by  the  column  are  greater  than  the  Euler  load 
when  the  loading  is  rapid,  the  column  is  slender  and  the  ini¬ 
tial  deflections  are  small.  Sevln  /4/  confirmed  the  previ¬ 
ous  results  while  including  the  effects  of  axial  inertia 
which  were  not  considered  by  Hoff.  He  Inferred  that,  for  an 
elastic  column,  axial  inertia  effects  are  of  negligible  im¬ 
portance  with  regard  to  the  gross  behavior  of  conventional 
structural  columns,  regardless  of  the  initial  deflected  shape, 
end  constraint,  or  type  of  axial  loading. 

Taylor  /5/  investigated  the  case  of  an  eccentrically 
loaded  elastic  column  subjected  to  a  half-sine  force  pulse 
and  he  included  the  effects  of  axial  inertia.  Using  a  limit¬ 
ing  value  of  extreme  fiber  strain  as  the  failure  criterion 
he  determined  the  column  failure  load  quantitatively.  He 
demonstrated  that  elastic  columns  will  support  loads  much 
greater  than  the  Euler  load  when  subjected  to  rapid  loading 
and  that  axial  inertia  effects  become  significant  with  this 
type  of  loading  as  the  rapidity  of  loading  is  increased. 

Gerard  and  Becker  /6/  studied  column  behavior  under  con¬ 
ditions  of  impact,  by  utilizing  the  properties  of  a  propa¬ 
gated  elastic  compressive  stress  wave  produced  on  failure 
of  a  tension  specimen.  They  concluded  that  a  column  can 
momentarily  support  a  dynamic  compressive  stress  of  any 
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magnitude  and  that  buckling  may  be  confined  to  a  length  less 
than  the  entire  length  of  the  column. 

The  inelastic  column  was  investigated  by  Chawla  /7/ 
in  essentially  the  same  manner  as  Hoff  (with  an  Initial 
curvature  or  effective  eccentricity)  and  he  found  that  for 
relatively  slow  rates  of  end  displacement  the  maximum  column 
load  was  less  than  the  tangent  modulus  load.  Brooks  and 
Vv’ilder  /8/  studied  the  effects  of  dynamic  loading  on  an 
Idealized  H-section  inelastic  column  with  an  initial  lateral 
deflection  using  a  constant  velocity  end  displacement.  From 
their  results  they  concluded  that  the  static  maximum  load 
may  be  employed  as  a  conservative  estimate  of  the  maximum 
load  of  a  column  regardless  of  the  rate  of  end  displacement. 

Investigations  of  the  dynamic  behavior  of  inelastic 
columns  are  limited  in  number  and  the  problem  of  determin¬ 
ing  the  behavior  of  inelastic  columns  subjected  to  pulse 
loading  appears  to  have  been  neglected.  Yvhereas  constant 
velocity  loading  leads  to  the  ultimate  failure  of  the  column, 
pulse  loading  may  not  necessarily  result  in  failure.  Pulse 
loading  allov.’3  the  column  to  recover  from  a  deflected  posi¬ 
tion  with  an  amount  of  residual  deflection  while  retaining 
a  capacity  to  sustain  subsequent  load.  The  problem  con¬ 
sidered  in  this  investigation  is  that  of  determining  the 
dynamic  behavior  of  a  pin-ended  inelastic  column  subjected 
to  a  half-sine  displacement  pulse  at  one  end  applied  with 
small  eccentricity.  The  problem  involves  the  prediction 


3 


of  the  maximum  force  developed  In  the  column  for  a  specified 
pulse  neight  and  duration  of  loading.  The  problem  also  In¬ 
cludes  the  prediction  of  the  residual  lateral  deflection 
resulting  from  the  pulse  loading.  The  effects  of  pulse  height 
and  duration  are  investigated  for  two  values  of  eccentricity 
and  three  column  slenderness  ratios. 
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2 •  DEVELOPMENT  OF_  THE  EQUATIONS  FOR  THE  MATHEMATICAL 

MO  DEI. 

The  column  considered  in  the  investigation  is  initially 
straight,  pin-ended,  and  of  constant  cross-section,,  The  ef¬ 
fects  of  axial  inertia  and  shear  strains  are  neglected,  but 
non-linear  axial  strain  components  due  to  bending  at  the 
center  portion  of  the  column  length  are  considered.  The 
effect  of  a  non-linear  stress-strain  curve  on  loading  and 
departure  from  this  curve  on  unloading  is  also  included. 

The  loading  on  the  column  is  produced  by  axial  displace¬ 
ment  of  the  pin  at  one  end.  The  loaded  end  is  free  to  rotate 
about  the  pin,  which  is  eccentrically  placed.  The  other  end 
is  constrained  by  a  fixed  pin  of  equal  eccentricity.  The 
end  displacement  is  imposed  by  a  moving  ram  executing  a  half¬ 
sine  displacement  pulse.  The  ram  can  exert  only  compressive 
force  on  the  column  (see  Fig  .  1). 

The  problem  is  developed  and  solved,  not  in  terms  of  the 
real  column,  but  in  terms  of  the  lumx>ed  parameter  mathematical 
model  shown  in  Fig.  1.  The  equations  for  the  model  are  then 
solved  in  a  high  speed  digital  computer.  The  force-deflection 
characteristic  of  the  springs  is  derived  from  the  material 
stress-strain  diagram.  The  assumed  behavior  is  indicated  in 
Fig.  2  where  the  solid  curve  OABD  represents  the  variation 
of  compressive  stress  (T  under  montonically  increasing  com¬ 
pressive  strain  6  .  If,  upon  reaching  point  B,  the  strain- 
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Pig.  1*  Geometry  of  the  Mathematical  Column  Model, 


I 
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ir.g  is  reversed,  the  path  will  be  along  BC  whose  slope  is 
that  of  the  tangent  to  the  solid  curve  at  0.  A  second  re¬ 
versal  will  cause  the  compressive  stress  to  increase  again 
following  CB.  After  again  reaching  the  curve  at  B,  the  path 
will  be  along  the  solid  curve  toward  D0  Additional  strain 
reversals  would  result  in  additional  "detours"  parallel 
to  BC.  In  the  study  the  stress-strain  relationship  after 
departure  from  the  compressive  stress- strain  curve  is  con¬ 
sidered  linear  into  the  tensile  stress  region.  The  maximum 
tensile  stresses  in  the  column  model  are  limited  to  values 
near  the  nominal  yield  stress  of  the  material.  It  is  pre¬ 
sumed  that  any  tensile  stress  in  the  column  model  springs 
is  insufficient  to  cause  plastic  straining. 

This  behavior  is  incorporated  in  the  springs  by  taking 
the  stiffness  of  each  as  K=A,  ET/JL ,  where  A,  is  half  the  col¬ 
umn  cross-sectional  area,  ET  is  the  tangent  modulus,  and 
is  the  unstrained  length  of  a  spring.  Again  referring  to 
Fig.  2,  at  any  point  on  the  solid  curve  Er  is  the  slope  of 
the  curve.  For  points  below  the  curve  E~=E  (Young's  Modulus). 

This  model  having  only  two  degrees  of  freedom  was  cnosen 
for  the  investigation  in  order  to  limit  the  computer  solution 
time  to  reasonable  values,  but  yet  to  retain  the  effects  of 
Inelastic  column  behavior.  It  is  felt  that  this  model  is  a 
good  compromise  in  that  it  approaches  the  deflected  shape 
of  the  inelastic  column  and  also  retains  the  effects  of 
lateral  inertia  and  of  the  longitudinal  strains  on  each  side 
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of  the  column  due  to  both  the  axial  load  and  bending. 

From  the  free  body  diagram  of  half  the  column  model 
shown  in  Fig.  3,  using  d'Alembert’s  principle,  the  follow¬ 
ing  equations  may  be  written: 


P  =  f;  +  Fz 


(1) 


Q—11 


dt* 


(2) 


P  (e  cos-©-  +  sine-)-  Q(±*.  cos©-  -  e  sin-e)  = 

(Fi cos^ 

Dividing  the  latter  by  cose  and  substituting  for  Q  from 
equation  (2) 


P(e 


ytan-0-)  - 


dt* 


—  e  "tan-e-)  - 


(3) 
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From  kinematic  considerations  of  the  column  model  in 
a  deflected  position  the  following  equations  may  be  written 


c£,  =  cC  +  2L^sin-e 


(4) 


dz=  c^m  -  z  sin© 


(5) 


where  c£,  and  are  the  axial  shortening  of  Spring  1  and 
Spring  2  respectively  and  c(m  is  the  average  axial  short¬ 
ening  at  the  center  of  the  model.  Further 


=  +  2.  •  4f-(j  -cos©)  +  z  e  sin© 


(6) 


y  =-jr  sin©  -  e(i-cos-e-) 


(7) 


where  x  is  the  end  displacement  at  the  loaded  end  of  the 
model  and  y  is  the  lateral  deflection  at  midlength  in  the 
model. 

The  following  approximations  are  consequences  of  assum 
ing  that  is  small  and  e/L  is  small: 

~  »  e  t  a  n  © 

sin-©  »  cos©) 

t  an-©  ~  sin©  ^  © 
l  -  cos  ©  »  4r  ~&Z 

An  evaluation  of  the  errors  resulting  from  these  approxima- 
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tions  is  given  in  App.  VI. 

Equation  (7)  now  becomes: 

7  =  (8) 

Now,  substituting  the  approximations,  together  with  equa¬ 
tion  (8),  into  equations  (3)  through  (6)  and  rearranging, 
results  in 

^  =  4ECe+y)-(F;-F,)^  (9) 

c£i  ■  dm  +  y 

d -  dm  -  y 

Solving  for  and  substituting  in  the  above  expressions 
for  c£,  and  gives 

<£,  =  x-  -  -try1  -  y  +  ^  y 
-ty1  -  ^-y  -*s§y 

Since  the  model  has  but  two  degrees  of  freedom,  only  a 
limited  correspondence  may  be  established  between  column  and 
model  parameters.  Energy  considerations  are  used  for  this 
purpose. 


(10) 

(11) 
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To  determine  this  correspondence  both  the  column  and  the 
model  are  considered  for  the  elastic  case.  Figure  4  indicates 


an  elastic  column  in  a  deflected  position  during  free  vibra¬ 
tion.  It  has  an  assumed  half-sine  deflected  shape  and  the 
motion  is  considered  sinusoidal  with  time. 

The  maximum  potential  energy  is  found  from  the  elastic 
strain  energy  due  to  bending  and  is 

,2. 


El 


Potential  Energy  = 


dx 


where  412  is  the  curvature  of  the  axis  of  the  column  and  I 
d  x2- 

is  the  moment  of  inertia  of  the  column  cross-sectional  area. 
I  is  equal  to  ^"A,  where  A  is  the  column  cross-sectional 
area.  This  can  be  shown  to  be 


Potential 


Energy  =  31 


4  El 
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(column) 


The  maximum  kinetic  energy  of  the  column,  neglecting 
axial  inertia  effects,  is 

Kinetic  Energy  =-|  C 

'o 

which  can  be  shown  to  be 


Kinetic  Energy 


=i 


A  (On  y 


(column) 


where  ^  is  the  mass  per  unit  length,  <jjn  is  the  circular 
frequency  of  first  mode  lateral  vibration  of  the  column  and 
Lis  the  column  length.  Setting  the  maximum  potential  energy 
equal  to  the  maximum  kinetic  energy  and  solving  for  the  circular 
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Fig.  4.  Elastic  Column  In  Free  Vibration. 
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frequency  of  first  mode  lateral  vibration  of  the  column  re¬ 
sults  in 


( column) 


Figure  5  shows  the  model  with  non-eccentric  pin  connec¬ 
tions  in  a  deflected  position.  The  stiffness  of  each  spring 
1 s  and  the  motion  of  the  mass  is  again  considered  sinusoi¬ 

dal  with  time.  The  maximum  potential  energy  as  found  from 
the  elastic  strain  energy  in  the  springs  is 


Potential  Energy  =  %  +  <£z  ) 


Considering  the  model  to  be  in  free  vibration  without  axial 
load  the  average  spring  deflection  is  zero  or  c£,  =  -  <^z  • 
Under  these  conditions  it  can  be  shown  from  equations  (10) 
and  (11)  that 


and  the  model  maximum  potential  energy  becomes 
Potential  Energy  =  1  .  (N  ^ — y* 


/m 


(model) 


r  \i 

The  maximum  kinetic  energy  of  the  model  neglecting  axial 
inertia  effects  is 


Kinetic  Energy  = 


(model) 


Now,  setting  the  maximum  potential  energy  equal  to  the  max¬ 
imum  kinetic  energy  and  solving  for  the  circular  frequency 
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of  first  mode  lateral  vibration  of  the  model  results  in 


(model ) 


To  determine  the  undeflected  length  of  the  model 
springs  the  potential  energies  of  the  column  and  model  are 
equated.  Taking  the  center  deflections  of  column  and 
model  to  be  equal,  Si.  is  found  to  be 


(12) 


To  determine  the  relation  between  column  and  model 
masses  the  kinetic  energies  of  the  column  and  model  are 
equated.  Again  taking  the  center  deflections  to  be  equal, 
m  is  found  to  be 


With  the  two  previous  correspondences  it  can  be  shown  that 
the  circular  frequencies  of  the  column  and  the  model  are  the 
same. 

Now,  considering  the  differential  equation  of  motion 
(as  previously  derived  for  the  model)  for  the  elastic  case, 
together  with  the  results  of  the  elastic  energy  development, 
Euler  load  and  end  displacement  correspondences  between  the 
column  and  model  are  determined.  These  will  be  used  later  in 
the  analysis  to  define  non-dimensi oral  load  and  displacement 
parameters.  If  the  acceleration  and  the  eccentricity  of  the 
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column  model  are  both  zero,  equation  (9)  becomes 


Py  =  (F,-F4)^ 


(») 


Equation  (14)  applies  to  static  buckling  of  an  elastic  column 
with  zero  eccentricity.  Now,  for  the  elastic  case 

and 

From  equations  (10)  and  (11) 

Thus 


Substituting  this  into  equation  (14),  the  Euler  load  Pe  of 
the  column  model  is  found  to  be 


P  =  T T4  A,  C 
e  4  (L/^)2- 

where  L/^  is  the  slenderness  ratio  of  the  column.  In  terms 
of  the  total  cross-sectional  area  the  equation  becomes 


p  _  JI4  A  Z 

c  ‘  s'  TL/O1 


(15) 


The  Euler  load,  in  terms  of  the  Euler  displacement,  is 

Pz~  ZK  Xe 


and,  substituting  for  K  and  PE 


the  Euler  displacement  becomes 

(16) 
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With  these  results  determined  from  elastic  behavior  of 
the  column  and  model,  the  equations  for  the  strains  and  for 
dynamic  equilibrium  of  the  inelastic  column  model  are  develop¬ 
ed  in  a  convenient  non-dimensional  parametric  form.  From 
equation  (10)  the  strain  in  Spring  1  is 


=  _cL-  —  nl  cL 

X  32.  L 


or 


f  s  -Hi 

3  2.  LX 


L3  J  1}  J  LT 

Defining  a  non-dimensional  lateral  deflection  as  V=  y/e  and 
making  the  substitution  y=Ve,  together  with  the  substitution 
for  the  Euler  displacement  from  equation  (16),  the  strain  may 
be  expressed  as 
■4 


r  —  — — — —  jl  _  (-£_fv  +  -fiL  V 

^  i  8  |_  Z  V  ^  J  v  V  ^  x  +  ^  J 


(17) 


Applying  this  same  reasoning  for  the  strain  in  Spring  2  re¬ 
sults  in 


TT 


J 


(18) 


From  these  equations  for  the  strains  at  the  center  por¬ 
tion  of  the  column  model  it  can  be  recognized  that  the  first 
terms  are  the  dimensionless  end  displacements,  the  second 
terms  are  the  corrections  to  this  end  displacement  due  to 
the  angular  rotation  of  the  rigid  rods,  the  third  terms  are 
the  corrections  due  to  the  change  of  the  effective  eccentri¬ 
city  with  lateral  deflection,  and  the  fourth  terms  result 
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from  bending 


The  relation  between  the  strains  and  stresses  during 
initial  loading  of  the  column  model  is  given  by  the  Ramberg- 
Osgood  /9/  equation: 


where  6  is  the  compressive  strain.  O’  is  the  compressive 
stress,  CT\j  is  the  yield  stress  of  the  material  and  n  is 
the  shape  parameter  of  the  stress-strain  curve.  By  defining 


a  non-dimensional  stress  S  as  the  previous  equation  may 

c. 

be  written  as 


€=  S 


(19) 


This  is  the  equation  which  is  used  to  determine  the  stress 
in  each  spring  as  the  model  is  initially  loaded.  For  the  com¬ 
pressive  stress-strain  relation  considered  in  this  investi¬ 
gation  the  shape  parameter  is  n^lO  and  the  dimensionless 
yield  stress  is  Sip0.004.  The  stress-strain  curve  is  shown 
in  Fig.  6  and  is  typical  of  24S-T  aluminum  alloy  material. 

As  each  of  the  springs  unloads,  it  leaves  the  stress- 
strain  curve  following  a  straight  line  of  slope  1  . 

As  subsequent  loading  takes  place  in  either  spring  the  same 
straight  line  is  retraced  and,  as  the  stress  reaches  a  level 
that  puts  it  on  the  stress-strain  curve  again,  it  follows 
the  curve  until  further  unloading  occurs. 

The  differential  equation  of  motion  is  developed  in 
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nondimenslonal  form  by  dividing  equation  (9)  by  the  eccentri¬ 
city  e  and  the  square  of  the  circular  frequency  of  first  mode 
lateral  vibration  COn  and  substituting  from  equation  (1)  for 
P.  After-  rearranging  terms  the  equation  becomes 


cT  C//e) 
^  d  t*- 


A  dimensionless  time  T  is  defined  as  T  =  (_Ont,  so  that 

(jJz  (dt)z  =  (dT)z  .  Then  — U-  becomes  ^  . 

d  tl  d 

By  making  the  substitution  V  for  y/e  this  becomes  which 

d  T 

♦  *  *  * 

is  symbolized  as  V.  By  substituting  V  and  V  and  recognizing 

p 

the  form  as  the  previously  defined  non-dimensional  stress, 

the  differential  equation  of  motion  becomes 


v = ^  (L/tf  [( s, + so(y+  o  -  (s,  -  st) 


(20) 


The  two  terms  in  the  brackets  can  be  recognized,  respective¬ 
ly,  as  proportional  to  the  external  bending  moment  due  to  the 
force  on  the  end  of  the  column,  and  the  internal  bending  mom¬ 
ent  due  to  the  spring  forces. 

In  order  to  solve  this  equation  it  is  necessary  to  specify 
the  forcing  function.  This  is  accomplished  by  specifying  the 
displacement  of  the  ram  (as  shown  in  Fig.  1 )  as  a  function  of 
time.  When  the  ram  is  in  contact  with  the  column  model  the 
ram  displacement  is  the  end  displacement  of  the  column.  If 
the  ram  is  not  in  contact  with  the  end  of  the  column,  the 
force  is  zero  (S|=-S^)  and  the  column  exhibits  free  vibration. 
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In  either  case  the  conditions  of  compatibility  and  dynamic 
equilibrium  are  maintained. 

The  ram  displacement  is  a  half-eine  pulse  of  height  X 
and  duration  .  To  enable  a  selection  of  dimensionless 
parameters  having  a  relationship  to  the  column  the  follow¬ 


ing  dimensionless  quantities  are  defined: 


(a) 


the  height  of  the  half-sine  ram  displacement  pulse 
as  a  fraction  of  the  Euler  end  displacement; 


(b)  the  ratio  of  the  natural  period  of  first  mode 
lateral  vibration  of  the  column  to  the  duration 
of  the  half-sine  pulse; 


With  these  dimensionless  quantities  and  the  previously  defined 


dimensionless  time.  The  dimensionless  ram  displacement 
becomes 


~X  ~  s  \  n  <y>T 


and 


for 


O  —  T 


^  2JX 
“  $ 


(21) 


"X.  =  O  for  T  > 

V 

To  determine  the  maximum  force  and  residual  lateral 
deflection  resulting  from  a  specified  ram  displacement  it  is 
necessary  to  compute  the  dynamic  response  of  the  model  to  the 
loading.  This  is  accomplished  through  a  numerical  integration 
process  applied  to  the  differential  equation  of  motion  and  a 
numerical  iteration  scheme  to  satisfy  the  relations  among  the 
strains,  stresses,  end  displacement  and  lateral  deflection. 
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3.  THE  NUMERICAL  SOLUTION 


The  equations  for  the  model  are  solved  utilizing  FOR¬ 
TRAN  programming  and  a  Control  Data  Corporation  1604  high 
speed  digital  computer.  The  investigation  is  accomplished 
with  two  basic  programs,  COL  1  and  COL  2.  The  former  is  em¬ 
ployed  primarily  to  obtain  graphical  output  directly  to  study 
the  dynamic  response  of  the  model  to  a  specified  loading.  Pro¬ 
gram  COL  2  is  employed  to  determine  the  numerical  values  of 
maximum  force  and  residual  lateral  deflection  for  a  specified 
loading  on  the  model.  The  complete  programs  in  FORTRAN  lan¬ 
guage,  together  with  the  program  notation,  are  given  in  Appen¬ 
dix  IV. 

The  numerical  integration  procedure  is  the  same  for 
both  programs,  but  the  program,  COL  2,  employed  to  determine 
maximum  values  is  more  comprehensive.  For  this  reason  the  de¬ 
tailed  discussion  is  limited  to  this  program.  A  discussion 
of  the  differences  between  these  programs  follows  the  detail¬ 
ed  discussion.  Figure  7  is  a  simplified  block  diagram  of 
Program  COL  2  and  Fig.  8  is  a  simplified  block  diagram  of  the 
subroutine  used  in  both  programs  to  determine  the  stress  from 
the  strain. 

Referring  to  Fig.  7,  the  input  to  the  program  includes 
the  following  dimensionless  parameters: 

(a)  L/^  (slenderness  ratio)  and  e/o  (eccentricity  ratio) 
to  prescribe  the  column  model. 

(b)  cx.  (ram  pulse  height)  and  (ratio  of  natural  period 
of  first  mode  lateral  vibration  to  ram  pulse  duration) 
to  prescribe  the  forcing  function  on  the  model. 
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Fig.  7*  Simplified  Block  Diagram  of  Computer  Program  COL  2 
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(c)  n  ( Ramberg-Osgood  shape  parameter)  and  Sy  (yield 
stress  of  the  material)  to  prescribe  the  compressive 
3tress-strain  curve, 

(d)  AT  (time  Increment  )  to  prescribe  the  Incremental 
values  In  the  numerical  Integration  process. 

The  Initial  conditions  are  then  set  to  start  the  numerical 

Integration  process. 

Henceforth,  the  dimensional  quantities  and  the  dimension¬ 
less  quantities  are  referred  to  synonomously  except  where  con¬ 
fusion  may  result,  i.e.,  lateral  deflection  for  dimensionless 
lateral  deflection,  etc. 

The  numerical  Integration  scheme  used  in  the  programs 
1 8  a  Runge-Kutta  process  /10/.  This  process  predicts  the 
lateral  deflection  V  and  the  lateral  velocity  V  at  a  time 
T+AT  from  four  estimates  of  the  lateral  acceleration  V  cal¬ 
culated  in  the  time  increment  AT,  Accuracy  Is  of  the  order 
of  (  A  T)3  .  Within  this  predictor  loop  further  calculations 
are  made  to  satisfy  the  relations  among  the  strains,  stresses, 
end  displacement  and  lateral  deflection. 

The  ram  displacement  Is  computed  from  equation  (21). 

The  values  of  ram  displacement  and  of  lateral  deflection  are 
then  used  to  determine  the  strain  in  Spring  1  from  equation 
(17).  The  main  program  then  calls  Subroutine  STRESS,  Because 
this  routine  is  written  with  dummy  variable  Inputs  to  accommo¬ 
date  the  strains  and  stresses  In  both  Spring  1  and  Spring  2, 
the  following  discussion  is  presented  based  on  a  general 
strain  €  and  a  general  dimensionless  stress  S. 


% 
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Subroutine  STRESS,  shown  in  Fig.  8,  is  entered  with  the 
present  value  of  strain  £  and  the  previous  values  of  stress  S 
and  strain  £a  calculated  for  the  preceding  time  T-AT.  The 
solid  curve  in  Fig.  9  represents  the  dimensionless  compres¬ 
sive  stress-strain  curve.  The  dashed  vertical  lines  indicat¬ 
ed  by  A,  B,  and  C  in  the  figure  are  representative  of  con¬ 
ditions  which  must  be  considered  in  arriving  at  a  value  for 
the  stress  S  from  a  value  of  the  strain  €  . 

Referring  now  to  Figs.  8  and  9,  if  the  present  value  of 
strain  is  negative  this  condition  is  represented  by  A.  The 
spring  has  unloaded  from  the  stress- strain  curve  following 
a  straight  line  with  a  slope  equal  to  unity  and  the  stress 
is  calculated  from 

S  =  $a  +  (,£  “  €a)  (22) 

This  equation  is  also  used  to  determine  the  stress  S  if  €.=0. 

If  £  is  positive  there  are  three  possibilities  to  be  consider 
ed„  These  are: 

(a)  loading  is  occurring  along  the  compressive  stress- 
strain  curve; 

(b)  unloading  is  occurring  on  a  unity  slope  line; 

(c)  loading  is  occurring  on  a  unity  slope  line. 

To  determine  which  of  these  possibilities  exists  the  incre¬ 
mental  change  in  strain,  A£=  6  -£a  ,  Is  calculated.  If 

this  quantity  is  negative  this  condition  is  represented  by 
B.  The  spring  is  unloading  on  a  unity  slope  line  and  the 
stress  is  calculated  from  equation  (22).  If  the  incremental 
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Fig,  8.  Simplified  Block  Diagram  of  Computer  Program  for 
Subroutine  STRESS* 


Plfi*  9*  Stress-Strain  Considerations  In  Computer  Program 
for  Subroutine  STRESS* 


/ 
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change  in  ©train  is  positive  two  possibilities  exist*  Either 
the  stress  is  on  the  compressive  stress-strain  curve  or  the 
stress  is  increasing  on  a  unity  slope  line*  This  situation 
is  represented  by  C  in  Fig*  9*  The  stress  level  is  deter¬ 
mined  by  calculating  a  value  of  stress  from  equation  (22) 
and  by  calculating  a  value  of  stress  from  the  Ramberg-Osgood 
equation  (19)  and  then  selecting  the  smaller  of  these  two 
values  as  the  stress  existing  in  the  spring*  If  the  incre¬ 
mental  change  in  strain  is  zero,  the  stress  S  is  equal  to  S3  , 

To  calculate  the  value  of  stress  from  the  Ramberg-Osgood 
equation  it  is  necessary  to  compute  the  real,  positive  root 
of  the  equation*  This  is  accomplished  by  a  Newton-Raphson 
Iteration  /10/  where  the  strain  f,  is  used  as  the  initial 
estimate  for  the  root  Sj.  *  The  stress  S  could  have  been  de¬ 
termined  from  the  strain  €  each  time  by  calculating  a  value 
from  equation  (22)  and  a  value  from  the  Ramberg-Osgood  equa¬ 
tion  (19)  and  selecting  the  smaller  value*  However,  this 
would  involve  an  iteration  for  each  stress  calculation,  thus 
increasing  the  computer  time*  Therefore,  this  method  was  not 
used* 

After  S  is  determined  the  value  for  the  strain  £z  is 
calculated  from  equation  (18)*  The  main  program  again  calls 
Subroutine  STRESS  to  determine  the  value  of  S*  „ 

At  this  point  the  program  calculates  the  sum  S,+Sz, 
which  is  directly  proportional  to  the  force  on  the  end  of 
the  model,  to  determine  if  it  is  compressive  (positive)  or 
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tensile  (negative).  As  previously  discussed  for  the  model 
configuration  (Fig.  1),  the  ram  cannot  transmit  a  tensile 
force  to  the  model.  Therefore,  if  this  sum  is  negative  the 
program  Initiates  an  iterative  Interpolation  scheme  to  deter¬ 
mine  the  end  displacement  of  the  model  which  will  cause  the 
sum  Si+Sj,  to  be  zero.  This  involves  recalculating  the  strains 
and  then  the  stresses  for  each  end  displacement  until  the  force 
is  zero.  The  model  is  then  in  a  state  of  free  vibration.  This 
free  vibration  continues  until  the  ram  again  contacts  the  column 
or  until  the  calculation  terminates. 

Once  these  states  of  stress  (end  force  either  positive 
or  zero)  are  determined,  the  lateral  acceleration  V  is  cal¬ 
culated  from  equation  (20).  After  four  passes  are  made  to 
this  point  in  the  program  to  determine  the  values  for  the 
Runge-Kutta  predictions,  as  previously  discussed,  the  lateral 
velocity  V  and  the  lateral  deflection  V  are  calculated 
for  the  end  of  the  time  interval  A  T. 

At  this  point  the  program  computes  the  force  which  can 
be  shown  to  be 


and  a  continuous  check  is  made  to  determine  the  maximum  value 
of  this  quantity  which  is  stored  for  print-out.  This  contin¬ 
uous  check  is  also  made  to  determine  the  maximum  value  of 
compressive  stress  in  each  of  the  springs.  These  values  are 
stored  for  use  in  strain  energy  calculations  to  be  made  later. 
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A  continuous  check  Is  also  made  to  determine  the  extreme  stress 
(tensile  or  compressive)  In  each  spring  for  prlnt-out.  The  pro 
gram  then  calculates  the  ram  work  done  on  the  model  for  the  pre 
ceding  time  Increment  and  adds  this  quantity  to  the  previous 
work  done  by  the  ram. 

This  whole  process  continues  until  the  ram  displacement 
pulse  has  terminated  and  the  column  Is  left  In  a  state  of 
free  vibration.  Then  a  continuous  check  Is  Initiated  to 
determine  the  time  at  which  S|=S2.=0.  The  lateral  deflec¬ 
tion  at  this  point  Is  the  residual  deflection  y>-  In  the 
model,  as  the  lateral  acceleration,  the  longitudinal  end 
force  and  the  Internal  moment  are  all  zero,  and  the  residual 
strains  determine  the  permanent  deformation  that  has  resulted. 

After  the  residual  deflection  has  been  calculated  a 
switch  is  set  to  detour  this  portion  of  the  program.  The 
new  path  Includes  a  continuous  check  to  determine  the  time 
at  which  the  lateral  velocity  is  zero  (point  of  extreme  lateral 
deflection  and  zero  kinetic  energy).  At  this  point  the  net 
energy  Input  to  each  of  the  springs  is  calculated  from  the 
area  under  the  stress-strain  curve  using  the  present  values 
of  stress  and  the  previously  stored  maximum  values  of  the  com¬ 
pressive  stress.  This  energy  is  a  combination  of  the  strain 
energy  of  plastic  deformation  and  of  the  potential  energy  due 
to  existing  elastic  strain  in  the  model.  This  energy  should 
be  equal  to  the  ram  work  done  on  the  model. 

When  the  energy  has  been  calculated  the  program  then 
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prints  all  input  quantities  and  the  results.  The  results 
include  the  values  of  P^/Pj.  ,y*-/^  ,Slrnax  ,  Sima%  and  the 
time  T  that  each  of  these  occur,  as  v.'ell  as  the  work  and 
energy.  The  times  are  included  to  complement  the  time  plots 
from  the  other  program  as  an  aid  in  visualizing  some  of  the 
dynamic  behavior  and  to  determine  the  problem  time  termina¬ 
tion  of  the  program. 

The  program  employed  to  obtain  graphical  output  directly 
is  Program  COL  1  (the  FORTRAN  language  listing  is  given  in 
Appendix  IV) .  This  program  differs  basically  from  Program 
COL  2  described  in  the  preceding  paragraphs  in  that  it  stores 
the  variable  quantities  for  each  time  T  in  array  form  instead 
of  storing  only  the  maximum  values.  These  array  quantities 
are  then  available  for  Subroutine  GRAPH  (not  shown  in  the 
listing  but  available  on  the  FORTRAN  Library  Tape  at  the  U.S. 
Naval  Postgraduate  School,  and  also  available  from  the  Control 
Data  Corporation).  Subroutine  GRAPH  is  designed  to  accept 
a  maximum  of  900  points  for  each  variable  to  be  plotted.  Thus 
the  program  is  written  to  calculate  only  900  points  while  main¬ 
taining  the  necessary  accuracy  for  the  numerical  integration 
process  compatible  with  the  time  increment  AT. 

In  addition  to  the  input  data  required  in  Program  COL  2, 
axis  scaling  and  multiple  curve  data,  graph  title  data  and 
curve  label  data  are  required  for  Frogram  COL  1.  In  addi¬ 
tion  to  the  graphical  output  the  program  prints  the  input 
quantities  and  the  values  of  T,  x/x£,  ,  S,  ,  Sz,  y/^ 

and  F/pEfor  every  fifteenth  point.  As  all  of  these  quantities 
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have  been  stored  as  9G0  point  arrays  any  combination  of  these 
variables  may  be  plotted  as  a  function  of  another.  Subroutine 
GRAPH  converts  these  900  point  arrays  to  input  data  on  a  mag¬ 
netic  tape  for  a  Control  Data  Corporation  16G  computer  which 
is  programmed  to  provide  output  as  input  to  a  plotter. 

An  example  of  one  set  of  curves  used  in  the  investi¬ 
gation  is  shown  in  Appendix  I.  The  example  includes  three 
graphs  with  two  curves  on  each  graph.  The  first  graph  is  a 
plot  of  €,  and  versus  T,  The  second  graph  is  a  plot  of 
St  and  versus  T  and  the  third  graph  is  a  plot  of  y/^  and 
F/Pe  versus  T,  Other  examples  of  graphical  output  are  shown 
in  Appendices  II  and  III.  However,  the  graphs  shown  in 
Appendix  II,  which  are  for  very  slow  loading  approaching  the 
static  case,  have  been  obtained  with  a  modified  version  of  the 
program.  The  modification  Includes  calculating  1800  points 
to  maintain  the  accuracy  for  the  numerical  integration  pro¬ 
cess,  then  re-storing  every  other  point  of  the  1800  point 
arrays  as  900  point  arrays  to  accomodate  Subroutine  GRAPH. 
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4 


RESULTS  AND  DISCUSSION 


a.  Character  of  the  response 

An  understanding  of  the  results  is  dependent  upon  an 
appreciation  of  the  distinctive  features  of  the  loading  and 
their  effects  upon  the  response.  Loading  Is  accomplished  by 
a  forced  displacement  of  one  end  of  the  column  in  the  form 
of  a  half-sine  pulse.  This  pulse  is  characterized  by  the 
dimensionless  parameters  cx.  and  „  The  quantity  is  the 
ratio  of  the  displacement  pulse  height  to  that  end  displace¬ 
ment  which,  in  the  absence  of  lateral  deflection,  corresponds 
to  the  Euler  buckling  load.  The  parameter  is  the  ratio  of 
the  first  mode  natural  period  of  lateral  vibration  of  the  col 
umn  model  to  the  pulse  duration.  Thus, ex  measures  the  ampli 
tude  of  the  loading  and  measures  its  rapidity. 

Consider  the  very  slow  loading  (Q}«l)  of  a  slender 
column  to  a  maximum  end  displacement  beyond  the  Euler  crlti- 
al  value  (oo  I  ).  At  axial  loads  above  the  Euler  limit 
the  straight  form  of  the  column  is  unstable.  To  initiate 
lateral  deflection  of  the  column  model  the  end  pins  are 
eccentrically  located.  The  effect  of  this  eccentricity  is 
to  produce  bending  at  loads  below  the  stability  limit  of  the 
column.  The  maximum  force  (Prn  )  developed  in  the  column 
becomes  less  than  the  Euler  load  (Pt  ).  However,  for  very 
slow  loading  the  maximum  force  developed  in  the  column  model 
is  the  static  maximum  load.  Increasing  values  of CX  above 
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that  corresponding  to  the  Euler  critical  value  do  not  change 
the  value  of  the  maximum  force  P^ . 

For  increasing  loading  rates  (increasing  values  of  ) 
the  lateral  motion  of  the  column  is  retarded  by  the  inertia 
of  the  mass  and  any  meaningful  correspondence  to  static  be¬ 
havior  is  lost.  However,  the  maximum  force  developed  in  the 
column  is  related  to  the  rate  of  loading.  Also,  under  pulse 
loading  the  inelastic  column  recovers  from  a  deflected  posi¬ 
tion  with  an  amount  of  residual  deflection  and  retains  a 
capacity  to  sustain  further  load.  This  capacity  may  be 
measured  in  terms  of  the  amount  of  residual  deflection 
which  results  from  the  loading.  This  residual  deflection  is 
equivalent  to  additional  eccentricity  for  any  subsequent  load¬ 
ing.  The  maximum  force  P^  and  the  residual  deflection  y r  appear 
to  be  the  most  important  performance  criteria  for  the  dynamic 
behavior  of  the  inelastic  column  subjected  to  a  pulse  loading. 
Dimensionless  forms  of  these  two  quantities  are  utilized  in 
the  graphical  form  of  the  results. 

b.  Range  of  variables 

Solutions  for  Pm/Pe.  (maximum  force)  and  for  yr/^  (residual 
deflection)  have  been  obtained  for  a  range  of  from  0.01  to 
1.0  and  for  a  range  of  cx.  dependent  on  the  column  slenderness 
ratio.  The  results  include  solutions  for  three  values  of  col¬ 
umn  slenderness  ratio  (L/^  )  and  two  values  of  eccentricity 
ratio  (e/^  ).  For  slenderness  ratios  of  70  and  50  the  mini¬ 
mum  pulse  height  considered  was  c<=  1.  For  a  slenderness 
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ratio  of  100  the  minimum  pulse  height  considered  was  cx=  2. 

For  values  of  cx less  than  these,  inelastic  action  is  insig¬ 
nificant,  The  maximum  values  of  oc.  considered  are  based  on 
the  maximum  value  of  tensile  stress  occurring  in  the  column 
model  springs.  For  the  investigation  the  stress-strain  re¬ 
lationship  after  departure  from  the  compressive  stress-strain 
curve  is  considered  linear  into  the  tensile  stress  region. 

For  this  reason  the  maximum  tensile  stresses  in  the  column 
are  limited  to  values  near  the  nominal  yield  stress  of  the 
column  material.  For  slenderness  ratios  of  70  and  50  the  max¬ 
imum  pulse  height  considered  was  cx.  =  2,5,  For  a  slenderness 
ratio  of  100  the  maximum  pulse  height  considered  was  cx=*  4. 

For  values  of  cx  greater  than  these  the  maximum  value  of  ten¬ 
sile  stress  occurring  in  the  column  model  is  considerably  in 
excess  of  the  nominal  yield  stress  for  some  rates  of  loading, 

c.  Maximum  force 

Figure  10  is  a  plot  of  families  of  curves  of  maximum  force 
Plri/PE  versus  for  a  constant  pulse  height  cx  .  The  families 
Include  slenderness  ratios  of  100,  70,  and  50  for  an  eccentri¬ 
city  ratio  of  0.05,  Figure  11  is  a  plot  of  the  same  families 
of  curves  for  an  eccentricity  ratio  of  0,01,  The  lover  limit 
of  ^  *  0,01  is  for  slow  loading  rates  approaching  the  static 
case.  The  upper  limit  of  Q  ~  1.0  is  approaching  the  region 
where  axial  inertia  effects  become  important  in  the  elastic 
column  as  reported  by  Taylor  /5/«  As  the  inelastic  column 
model  does  not  include  the  effects  of  axial  inertia  =  1.0 
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Fig.  11.  Values  of  maximum  force  Pm/PE  vs.  (  c<  constant,  e/ 
for  column  slenderness  ratios  L/<^  of  100,  70,  and  50 


1 8  the  upper  limit  considered  in  this  study. 

Figures  10  and  11  show  that  as  the  rate  of  loading  is 
increased  the  maximum  force  developed  in  the  column  increases. 
This  is  due  to  the  delay  in  the  development  of  lateral  deflec¬ 
tion  because  of  lateral  inertia.  These  figures  also  indicate 
that  as  the  rate  of  loading  is  decreased  the  maximum  force  de¬ 
veloped  in  the  column  approaches  a  limit  asymptotically.  In 
Fig,  10  this  limit  may  be  observed  for  each  slenderness  ratio 
as  all  of  the  curves  converge  to  a  single  value  at  Q  =  0.01. 
This  limit  is  the  static  load  capacity  of  the  column  for  an 
eccentricity  ratio  of  0.05.  In  Fig.  11  this  limit  is  not 
attained.  The  spread  of  points  indicates  that  for  this  eccen¬ 
tricity  (e/(>  =  0.01)  there  is  considerable  dynamic  effect  at 
Q  =  0.01.  However,  the  curves  of  minimum  cx.  f0r  each  slen¬ 
derness  ratio  at  =  0,01  appear  to  approach  a  limit  of  the 
static  load  capacity  of  the  column  fcr  an  eccentricity  ratio 
of  0.01.  A  comparison  of  these  figures  confirms  that  as  the 
eccentricity  increases  the  maximum  force  developed  in  the 
column  decreases. 

d.  Slow  loading 

It  is  possible  to  make  some  observations  concerning 
the  values  obtained  at  the  low  end  of  the  range  by  con¬ 
sidering  the  static  case  of  a  centrally  loaded  inelastic 
column.  For  this  ca3e  the  maximum  load  the  column  can  sup¬ 
port  is  between  the  tangent  modulus  load  PT  and  the  reduced 
modulus  load  P^  .  In  terms  of  the  Euler  load  the  dimensionless 
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tangent  modulus  load  Is 


(23) 


where  ET  is  the  tangent  modulus  or  the  local  slope  of  the 
compressive  stress-strain  diagram.  The  tangent  modulus  de¬ 
fines  the  value  of  stress  at  which  bending  begins  for  the 
static  case  of  a  centrally  loaded  Inelastic  column.  In 
terms  of  the  model  constants  developed  previously  this 
critical  stress  is 


0”cr  — 


_  TT4  E- 


8  (L_/^)4 

or  in  terms  of  dimensionless  stress  this  becomes 


Sr  -  -Hi  ( _ 1— ' \ 

0cr  ~  ft  T  \L/^i 


(24) 


The  ratio  of  Young's  modulus  to  the  tangent  modulus  in  terms 
of  the  Bamberg- Osgood  parameters  is 


X  =  l+nig*lL(^_) 


(25) 


where  n  =  10  (the  shape  parameter)  and  =  C.0C4  (the  dimen¬ 
sionless  compressive  yield  stress)  are  the  parameters  used 
in  the  investigation.  By  choosing  a  value  of  Scy-  and  sub¬ 
stituting  this  for  S  in  equation  (25)  the  ratio  E/ET  is  de¬ 
termined.  By  substituting  the  values  of  this  ratio  and  SCr 
in  equation  (24)  the  slenderness  ratio  is  determined.  In 
this  way  a  value  of  the  ratio  Ej/E  is  determined  correspond¬ 
ing  to  each  of  the  slenderness  ratios  100,  70,  and  50.  This 
value  of  Sr/S  equals  PT  /Pe  from  equation  (23). 
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In  terms  of  the  Euler  load  the  dimensionless  reduced 


modulus  load  is 


Pe  £ 


(26) 


where  ER  is  the  reduced  modulus.  The  reduced  modulus  de¬ 
fines  a  maximum  value  of  stress  which  can  be  approached  but 
not  exceeded  for  the  static  case  of  a  centrally  loaded  inel¬ 
astic  column.  In  terms  of  the  model  constants  developed  pre¬ 
viously  this  critical  stress  is 

rr  -  t4 
^  "  8  (L^ 

or  in  terms  of  dimensionless  stress  this  becomes 


S,r.  =  Oil 

°cr  8  E  Il/P 


la>; 


(27) 


The  reduced  modulus  13  dependent  on  the  tangent  modulus  and 
the  geometry  of  the  column.  For  the  column  model  it  can 
be  shown  that  the  reduced  modulus  is 
r  -  iLh 

t'K  "  E  +  Et 

or  In  dimensionless  form  this  becomes 


Ea.  zEt/e 

e  *  +  et/e 


(28) 


and  thl3  value  is  equal  to  PR/Pe  from  equation  (26).  By 
choosing  a  value  of  Sc>-  and  substituting  this  for  3  in  equa¬ 
tion  (25)  the  ratio  Et/E  Is  determined.  The  ratio  Et/E  also 
determines  E*/E  from  equation  (28).  By  substituting  the  values 
of  Er/S  and  3Cr  in  equation  (27)  the  slenderness  ratio  is  deter¬ 
mined.  In  this  way  a  value  of  the  ratio  ER/S  is  determined 
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corresponding  to  each  of  the  slenderness  ratios  100,  70,  and 
50.  This  value  of  E^/E  equals  P*/PE  from  equation  (26). 

The  values  of  the  tangent  modulus  load  and  the  reduced 
modulus  load  determined  from  these  two  methods  are: 

for  L/^  =  100,  PT/P£  =  1.000  and  PR/Pe  =  1.000; 

for  L/^  =  70,  PT/Pt  =  0.969  and  PR/PE  =  0.977; 

for  L/<^  =  50,  Pr/Pf-  =  0.649  and  PR/Pp_  =  O.69I  . 

From  Figs.  10  and  11  It  is  seen  that  In  most  cases  the  max¬ 
imum  values  at  Q  =  0.01  are  less  than  the  tangent  modulus 
load  although  the  spread  of  points  for  an  eccentricity  ratio 
of  0.01  indicates  that  dynamic  effects  are  more  pronounced 
the  smaller  the  eccentricity.  For  tills  reason  some  of  tne 
maximum  values  are  greater  than  the  reduced  modulus  load. 

To  observe  further  the  effect  of  eccentricity  on  the 
load  capacity  of  the  column  at  very  slow  loading  rates, a 
study  was  made  to  determine  tbs  behavior  of  the  column  under 
a  ramp  displacement  of  the  end  at  a  loading  rate  correspond¬ 
ing  to  that  available  to  commercial  testing  machines.  In 
terms  of  the  parameters  previously  used  to  describe  the  pulse 
loading  of  the  column  the  ramp  displacement  of  the  head  is 

X  =  Q>  T 

where  cx.  =  1.0  and  =  0.01.  Computer  solution  results  in¬ 
cluding  a  graphical  presentation  of  Strain  1  (  £,  )  and  Strain 
2  (  6^  )  versus  force  (P/Pe)  were  obtained  for  the  three  column 
slenderness  ratios  of  100,  70,  and  50  for  eccentricity  ratios 
of  C.05,  0.01,  and  0.001.  The  graphical  results  shown  in  App.II 
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exhibit  a  general  resemblance  to  the  experimental  strain 
versus  force  curves  of  actual  column  test  data  displayed  by 
Shanley  /l/.  From  the  results  given  in  Table  1  it  is  seen 
that  for  eccentricity  ratios  of  C.05  and  0.01  the  maximum 
load  is  less  than  the  tangent  modulus  load  for  all  three 
slenderness  ratios  considered.  However,  the  maximum  load 
increases  with  a  decrease  in  the  eccentricity  ratio.  As 
the  eccentricity  ratio  is  decreased  to  0.001  the  maximum 
load  is  greater  than  the  tangent  modulus  load  but  less  than 
the  reduced  modulus  load  for  a  slenderness  ratio  of  50.  For 
a  slenderness  ratio  of  70  there  is  considerable  dynamic  effect 
and  the  maximum  force  exceeds  the  reduced  modulus  load  by  a 


Table  I  Computer  and  Theoretical  Static  Load  Comparison 
(For  Ramp  Loading  of  Column  Model) 


L/<^  =  100  ' 

L/^>  =  70 

L/^>  =  50 

Static 

P-r/PE 

1.000 

0.969 

0.649 

Static 

PR/P£ 

1.000 

0.977 

0.691 

Dynamic 

( Computer) 

for 

e/^  =  0.05 

0.952 

0.839 

0.584 

for 

e/<^  =  0.01 

0.983 

0.916 

0.633 

for 

e/(>  =  0.001 

1.000 

0.981 

c  .665 
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small  amount.  For  a  slenderness  ratio  of  100  there  is  consider¬ 
able  dynamic  effect  also  but  the  solution  oscillates  about  an 
average  maximum  value  of  P^/  PE  =  1,0. 

The  tabulated  values  of  Pm/PE  for  an  eccentricity  ratio 
of  0.05  for  a  ramp  displacement  loading  are  in  agreement  with 
the  maximum  force  values  Pm/PE  for  pulse  loading  at  Q>  =  0.01 
in  Fig.  10.  The  strain  (  £,  and  £4)  versus  force  (P/PE )  curves 
for  the  ramp  displacement  loading  for  an  eccentricity  ratio 
of  0.05  shown  in  App.  II  display  no  dynamic  behavior.  In  Fig. 

11  (e/(D  =  0.01),  where  the  dynamic  effects  are  more  pronounced 

at  Q>  *  0.01  for  pulse  loading,the  forces  on  the  C*.=  l  curves 
for  slenderness  ratios  of  70  and  50  agree  with  these  tabulated 
results.  For  a  slenderness  ratio  of  100  the  smallest  value  of 
c<»  considered  for  pulse  loading  was  2  and  this  curve  appears 
to  approach  a  value  P^/P^  =  0.983  asymptotically.  The  strain 
versus  force  curves  for  the  ramp  displacement  loading  for  an 
eccentricity  ratio  of  0.01  shown  In  App.  II  display  an  insignif¬ 
icant  amount  of  dynamic  behavior  which  does  not  affect  the 
values  of  maximum  force  P_/P..  .  Therefore,  these  values  of 

P^./Pj.  may  be  considered  the  static  load  capacity  of  the  inel¬ 
astic  column  for  eccentricity  ratios  of  0.05  and  0.01. 

Appendix  III  contains  a  compilation  of  the  dynamic  re¬ 
sponse  curves  stress  (Si  and  3?  )  versus  time  (T)  for  a  range 
of  from  C.l  to  1.0  and  for  an  eccentricity  ratio  e/c^  =  0.01. 
This  graphical  compilation  includes  the  solutions  for  the  min¬ 
imum  and  maximum  values  of  pulse  height  consldex^ed  in  the 
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study  for  all  three  slenderness  ratios 0  Included  with  the 
compilation  Is  a  discussion  of  the  more  prominent  features 
of  the  curves,. 


e.  Residual  deflection 

Figures  12 ,  13,  and  14  are  famlles  of  curves  showing 


residual 


height  ex.  and  for  e/^  =  0o01  and  0o05o  Figure  12  Is  the 
family  for  L/C  =  100,  Figc  13  Is  the  family  for  L/^  =  70, 


and  Flg0  14  Is  the  family 


The  points  shown 


on  these  curves  are  from  the  identical  computer  solution 
results  as  the  points  shown  in  Figs0  10  and  11 » 


From 


shown  that  In  general 


the  values  of  residual  deflection  Increase  with  the  rapidity 
of  loading  to  a  maximum  value  then  decrease  abruptly „  This 
abrupt  decrease  is  because  of  the  delay  in  the  development  of 
lateral  deflection  due  to  lateral  Inert laQ  From  Fig.  13 
for  L/x^  =  70  and  Fig.  14  for  L/c^  =  50  it  is  shown  that  in  gen¬ 
eral  the  values  of  residual  deflection  are  essentially  con¬ 
stant  as  the  rapidity  of  loading  increases  until  an  abrupt 
decrease  occurs  due  to  the  effect  of  lateral  inert la0  The 
exception  for  cx  =  l  and  L/£  =  70  is  due  to  the  small  amount 

of  Inelastic  action  that  takes  place  for  this  pulse  height 
at  all  values  of  0}  „  The  amount  of  residual  deflection  may 
be  expressed  in  terms  of  the  amount  of  inelastic  action  that 
occurs.  From  the  graphical  compilation  of  stress  versus  time 


curves  in  Appc  Til  and  the  data  complied  in  App0  V  the  amount 


Fig.  12.  Values  of  residual  deflection  yy/^  vs. 

(cx.  constant,  e/^  °»0.01  and  0.05)  for  a 
column  slenderness  ratio  L/^  of  100 
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Fig.  13.  Values  of  residual  deflection  y*-/P  vs.  Q 
(cx.constant,  e/^*  0.01  and  0.05)  for  a 
column  slenderness  ratio  L/^  of  70. 
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Fig.  14.  Values  of  residual  deflection  yr/$  vs. ^ 
(  ck.  constant ,  e/<£=  0.01  and  0.05;  for  a 
column  slenderness  ratio  L/^  of  50. 


of  Inelastic  action  may  be  determined 0  For  a  slenderness 
ratio  of  100,  inelastic  action  in  the  column  is  limited  to 
maximum  compressive  stress  levels  below  the  value  of  the 
yield  stress  for  all  solutions  considered  in  the  investiga¬ 
tion.  For  a  slenderness  ratio  of  50,  inelastic  action,  in 
general,  occurs  at  much  higher  maximum  compressive  stress 
level s-near  or  above  the  yield  stress.  For  a  slenderness 
ratio  of  70,  inelastic  action  occurs  at  maximum  compressive 
stress  levels  Intermediate  to  those  for  slenderness  ratios 
of  100  and  50. 

A.  comparison  of  Fig.  12,  13,  and  14  Indicates  that  as 
the  slenderness  ratio  decreases  the  effect  of  eccentricity 
on  the  values  of  residual  deflection  is  less  pronounced. 
However,  decreasing  eccentricity  causes  a  shift  of  the  re¬ 
sidual  deflection  curves  to  longer  pulse  durations  (smaller 
values  of  §  ).  This  is  because  dynamic  behavior  is  more 

pronounced  the  smaller  the  eccentricity  and  lateral  Inertia 
effects  become  more  important  at  slower  loading  rates.  It 
may  also  be  observed  from  these  residual  deflection  curves 
that  the  maximum  values  of  residual  deflection  do  not  necessa¬ 
rily  occur  at  the  slowest  loading  rates.  The  minimum  values 
of  residual  deflection  generally  occur  at  the  fastest  loading 
rates.  The  only  exceptions  are  for  the  columns  of  largest 
slenderness  ratio  (L /o  =  100)  and  largest  eccentricity 
(e/^  =  0.05).  However,  the  residual  deflections  appear  to 
follow  the  same  trend  but  are  limited  by  the  maximum  value  of 
the  parameter^  =  1.0(  shortest  pulse  duration  considered  in  the 
investigation) . 
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5 


CONCLUSIONS 


From  the  results  discussed  in  the  preceding  section, 
the  following  conclusions  may  be  drawn  for  the  rates  of 
loading  considered  in  the  investigation: 

(a)  The  inelastic  column  model,  together  with  a  high 
speed  digital  computer  and  the  peripheral  plotting 
equipment,  provides  a  sound  and  reasonable  method 
for  the  study  of  the  dynamic  behavior  of  inelastic 
columns. 

(b)  For  very  slow  loading  (long  pulse  duration)  the  model 
exhibits  behavior  in  agreement  with  Shanley's  postu¬ 
lates. 

(c)  As  eccentricity  is  reduced,  the  lateral  response  is 
slowed.  Accordingly,  columns  with  very  small  eccen¬ 
tricity  show  significant  dynamic  effects  at  quite 
slow  loading  rates. 

(d)  Within  the  range  of  pulse  durations  studied,  the 
maximum  load  Increases  monotonically  with  decreas¬ 
ing  pulse  duration. 

(e)  The  ratio  of  maximum  dynamic  load  to  static  load 
capacity  decreases  as  the  slenderness  ratio  is 
reduced. 

(f)  Within  the  range  of  pulse  durations  studied,  the 
smallest  residual  deflections  generally  accompany 
the  shortest  pulse  durations. 
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APPENDIX  I 


Example  of  Graphical  Output  from  Program  COL  1 
Graphical  Notation 

The  computed  points  on  the  graphs  are  the  values  of  the 
dimensionless  variables  developed  In  the  basic  equations. 

The  graphical  notation  Is  as  follows: 


Symbol  on  Graph 

Corresponding  Symbol 
the  basic  equations 

SLENDERNESS 

L/<^ 

ECC 

cv 

<1) 

ALPHA 

C< 

BETA 

§ 

STRAIN  1 

C 

STRAIN  2 

STRESS  1 

s 

STRESS  2 

Si 

DEFLECTION 

y/^ 

FORCE 

p/pE 

TIME 

T 

Each  of  the  curves  on  the  graph  Is  labeled  to  correspond 
to  the  terms  In  the  graph  title.  The  x-axis  Is  the  axis 
of  abscissas  and  the  y-axis  Is  the  axis  of  ordinates.  The 
numbers  specifying  the  coordinates  are  integer  numbers  locat¬ 
ed  to  the  right  of  the  coordinate  on  the  x-axis  and  above  the 
coordinate  on  the  y-axis„  The  scale  quantity  Is  the  value  of 
the  first  marked  coordinate  from  the  origin  in  either  the 
x-direction  or  the  y -direct ion  expressed  in  powers  of  ten 
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to  locate  the  decimal  point .  The  symbol  E  In  the  scale  value 
indicates  that  the  number  following  is  the  exponential  power 
of  ten  to  be  applied  to  the  basic  scale  value,  i.e.,  1.00E+C1 
is  1.00  x  10+l  . 

Discussion  of  the  curves 

On  pages  1-5  to  1-7  following  this  discussion  are  shown 
curves  of  6, »  ,  S,  ,  Sz  ,  y/^  and  P/PE  versus  T.  These 

curves  are  the  solutions  for  the  inelastic  column  model  dy¬ 
namic  response  to  a  ram  pulse  loading  with  the  following  par¬ 
ameters  : 

(a)  slenderness  ratio,  L/^  =  100 

(b)  eccentricity  ratio,  e/^  =  0.05 

(c)  ram  pulse  height, cx.  =  4.0 

(d)  ram  pulse  duration  in  terms  of  =  C.2, 

T*  =  2  TT/Q>  =  10T T 

The  discussion  includes  some  of  the  more  prominent  features 
of  the  dynamic  behavior  of  the  column  model  under  ram  pulse 
loading. 

A  comparison  of  the  strain  versus  time  curves  with  the 
stress  versus  time  curves  indicates  the  amount  of  inelastic 
action  that  takes  place  in  Spring  1  and  Spring  2.  In  Spring  1 
the  compressive  stress  level  increases  above  that  in  Spring  2 
as  the  amount  of  bending  increases.  The  internal  bending 
moment  in  the  column  model  is  directly  proportional  to 
S,  -  Sg_  .  As  the  amount  of  bending  increases  Spring  2  de¬ 
parts  from  the  compressive  stress-strain  curve  along  a  straight 
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line  with  unit  elope  and  continues  into  the  tensile  stress 
region.  The  values  of  stress  in  Spring  2  remain  on  this 
straight  line  for  the  remainder  of  the  solution.  The  com¬ 
pressive  stress  level  in  Spring  1  continues  to  increase  to  a 
point  where  it  also  departs  from  the  compressive  stress-strain 
curve  along  a  straight  line  as  the  internal  bending  moment  in¬ 
creases  and  the  force  decreases  to  zero.  This  force  is  dir¬ 
ectly  proportional  to  S  +  S  0  When  the  force  becomes  zero 
the  end  of  the  column  has  departed  from  the  ram  which  exhibits 
a  half-sine  pulse  displacement.  The  column  model  is  then  in 
a  state  of  free  vibration  until  the  end  again  makes  contact 
with  the  ram.  The  compressive  stress  level  in  Spring  1  again 
returns  to  values  on  the  compressive  stress-strain  curve. 

This  alternate  departure  and  contact  between  the  column  and 
the  ram  continues  until  the  ram  displacement  pulse  nears  its 
termination.  The  column  then  continues  in  a  state  of  free 
vibration  where  S,  =  -  Si0  The  amount  of  residual  deflec¬ 
tion  may  then  be  determined  from  the  mean  value  of  the  deflec¬ 
tion  or  from  the  deflection  at  the  time  when  S  =  S-  =  0,  At 
this  point  the  internal  forces  and  the  Internal  bending  moment 
in  the  column  model  are  zero. 

The  effect  of  eccentricity  on  the  dynamic  behavior  of 
the  column  model  with  all  other  parameters  remaining  the 
same  may  be  observed  by  comparing  the  stress  versus  time 
curves  in  this  appendix  with  the  stress  versus  time  curves 
on  page  III-13  where  the  eccentricity  ratio  is  0,01,  It  may 
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be  noted  that  the  essential  differences  In  these  curves  are 
that  the  maximum  compressive  stress  levels  attained  In  Spring  1 
and  Spring  2  are  greater  for  the  smaller  eccentricity „  Also, 
contact  between  the  ram  and  the  column  terminates  sooner  for 
the  smaller  eccentricity 0 
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APPENDIX  II 


Theoretical _  and  Computer  Solution  Graphical  Comparisons  of 
Maximum  Loads  for  the  Static  Case 

The  notation  for  the  graphs  displayed  in  this  appendix 
is  explained  In  Appc  Ic 

The  graphB  on  pages  II-3  to  11=11  folowing  this  dis¬ 
cussion  are  the  curves  ?  and  •  ,  versus  P/P  for  the  three 
slenderness  ratios  considered  in  the  study ,  As  previously 
mentioned,  the  loading  for  the  study  in  this  appendix  is  a 
ramp  displacement  of  the  end  of  the  column  model  at  a  dis¬ 
placement  rate  comparable  to  that  in  commercial  testing  ma¬ 
chines,,  The  study  includes  eccentricity  ratios  of  0,05, 
0,01,  and  0o00l  for  each  of  the  slenderness  ratios,,  The 
tabulated  results  of  this  study  are  shown  in  Table  I  in 
Section  4, 

In  addition  to  the  previous  comments  made  concerning 
the  results,  it  may  also  be  seen  from  the  graphs  that  for 
all  slenderness  and  eccentricity  ratios  considered  an  amount 
of  bending  takes  place  before  strain  reversal  occurs  in 
Spring  2  and  before  the  maximum  load  is  reached.  This  was 
essentially  Shanley0s  point  in  proving  the  fallacy  of  the 
reduced-modulus  theory  for  column  action.  As  stated  by 
Shanley  /!/ 8 

,,,ocothe  reduced-modulus  theory  is  not  correct 
for  predicting  the  maximum  load  up  to  which  a  per¬ 
fect  column  will  remain  straight.  This  is  because 
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it  is  possible  for  the  column  to  bend  simultaneously 
with  increasing  axial  loada  Under  such  conditions 
it  is  possible  to  have  bending  without  introducing 
any  strain  reversal,  upon  which  the  reduced-modulus 
theory  depends 0 
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APPENDIX  III 


Column  Model  Dynamic  Response  Curves 

The  notation  for  the  graphs  displayed  in  this  appendix 
is  explained  in  App.  I. 

The  graphical  compilation  on  pages  III-3  to  III-26 
following  this  discussion  gives  the  cur  ves  St  and  Sz  versus 
T  for  Q  =  0.1,  0.2,  0.5  and  1.0  and  for  an  eccentricity  ratio 
e/^  =  0.01.  Included  in  the  compilation  are  the  solutions 

for  the  minimum  and  maximum  values  of  ram  pulse  height  cx.  con¬ 
sidered  in  the  investigation  for  all  three  slenderness  ratios. 
The  graphical  solutions  are  grouped  by  slenderness  ratio  with 
the  solution  for  the  maximum  value  of  cx.  followed  by  the  solu¬ 
tion  for  the  minimum  value  of  cx»  for  a  single  value.  The 
order  of  the  (^>  values  is  C.l,  0.2,  0.5  and  1.0,  that  is,  from 
the  longest  pulse  duration  of  loading  to  the  shortest  pulse 
duration  of  loading  considered  in  the  compilation. 

A  discussion  of  some  of  the  prominent  features  of  the 
dynamic  behavior  of  the  column  model  under  ram  pulse  loading 
is  included  in  App.  I  where  the  curves  of  6,  ,  6^,  S,  ,  Sz,  y/^ 
and  P./PE  versus  T  are  displayed  for  one  computer  solution. 

The  following  Is  a  discussion  of  some  of  the  more  prominent 
features  of  the  dynamic  behavior  of  the  column  model  in  terms 
of  the  stress  curve  compilation  in  this  appendix. 

In  general  the  amount  of  inelastic  action  in  the  column 
model  is  Inversely  related  to  the  slenderness  ratio.  For  a 
slenderness  ratio  of  100  inelastic  action  is  limited  to  cora- 
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pressive  stress  levels  below  the  value  of  the  yield  stress. 
For  a  slenderness  ratio  of  50  the  maximum  conpresslve  stress 
levels  are  near  or  above  the  yl9ld  stress.  For  a  slenderness 
ratio  of  70  inelastic  action  occurs  at  maximum  compressive 
stress  levels  intermediate  to  those  for  slenderness  ratios 
of  50  and  100. 

It  may  also  be  noted  that,  in  general,  contact  is  lost 
between  the  ram  and  the  end  of  the  column  (when  S(  =  -  S^) 
before  the  loading  pulse  nears  its  termination  for  values 
of  0.2  and  0.5  and  slenderness  ratios  of  70  and  50.  For  a 
slenderness  ratio  of  100  this  occurs  at  values  of  (^>  =  0.1, 

0 .2  and  0.5. 

It  may  be  seen  that  the  value  of  maximum  tensile  stress 
in  Spring  2  exceeded  the  value  of  the  nominal  yield  stress 
(3y  =  C.CC4)  for  the  values  of  0.2  and  0.5  and  ex  =»  2.5 
for  a  slenderness  ratio  of  50. 
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APPENDIX  _IV 

Comx'uter .Programs  and  Notation 
a.  General 

The  FORTRAN  language  computer  programs,  together  with 
the  notation,  and  the  program  operating  instructions  for 
both  Program  COL  1  and  Program  COL  2  are  included  in  the 
following  pages.  The  discussion  and  a  listing  of  Program 
COL  1  are  Included  on  pages  IV  -  2  to  IV  -10.  The  discuss¬ 
ion  and  a  listing  of  Program  COL  2  are  included  on  pages 
IV  -11  to  IV  -20.  The  FORTRAN  language  format  for  each  of 
the  programs  is  in  standard  FORTRAN  language  to  be  used  with 
a  standard  FORTRAN  compiling  routine.  The  program  language 
is  compatible  with  any  of  the  large  machine  FORTRAN  compiling 
systems  including  the  I3M  7090  and  the  CDC  1604  systems.  Each 
line  in  the  listing  is  toe  information  contained  on  one  80  col¬ 
umn  IBM  card  where  columns  1-5  contain  any  statement  number, 
column  6  is  used  to  indicate  a  continuation  card,  columns 
7-72  contain  the  FORTRAN  statement,  and  columns  73-80  are 
used  for  card  identification  only.  Although  the  notation  is 
basically  the  same,  the  subscript  notation  is  different  for 
each  program.  Program  COL  1  stores  a  number  of  variable 
quantities  in  array  form  for  print-out  and  to  accommodate 
Subroutine  GRAPH,  whereas.  Program  CCL  2  stores  only  one 
value  of  tne  variable  quantities  for  print-out. 
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be  Notation  for  Program  COL  1 

In  the  notation  for  Program  COL  1  the  subscript  J  re¬ 
fers  to  quantities  within  the  Runge-Kutta  prediction  loop 
which  are  calculated  for  intermediate  times  within  the  time 


Increment 

AT.  The  subscript  I 

refers  to  the 

quant  it les 

corresponding  to  each  time  T(I). 

The  notation 

for  Program 

COL  1  is 

as  follows: 

Computer  program  notation  for  Program 

C0L 1 

Computer 

Program 

Term 

Correspond ing 
terms  in  table 
of  symbols 

Computer 
Program 
term  _ 

Corresponding 
terms  in  tabl- 
of  symbol 3 

A 

n 

SI  GY 

S, 

ALPHA 

CX. 

SIG1( J) 

S, 

BETA 

SIG2( J) 

DDV ( J) 

V 

SLNESS 

L/^ 

DEFL(I) 

V 

STRNl(I) 

€, 

DELT 

AT 

ST RN2(  I) 

DISPL(I) 

x/xE 

ST  RSI ( I) 

3 1 

ECC 

eA> 

ST  RS2 ( I ) 

Sz 

EPSI1( J) 

T(r) 

r 

EPSI2( J) 

£-Z 

VSL 

V 

FORCE ( I ) 

P/PE 

X(J) 

x/xE 

o.  Operating  instructions  for  Program  COL  1 

The 

output  from  Program  COL 

1  includes  the 

print-out 

and  the  array  quantities  necessary  for  Subroutine  GRAPH. 
The  listing  for  Program  CCL  1  does  not  Include  Subroutine 
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GRAPH.  However,  a  listing  of  Subroutine  GRAPH  and  instruc¬ 
tions  for  its  use  may  be  obtained  from  the  Control  Data  Cor¬ 
poration.  To  use  the  program  to  obtain  the  print-out  quanti¬ 
ties  without  the  graphical  output  remove  the  cards  CALL  GRAPH. 
To  operate  Program  CGL  1  without  Subroutine  GRAPH  use  the  cards 
including  Subroutine  STRESS  in  the  order  shown  in  the  listing 
followed  by  one  blank  card,  then  the  input  data  card. 

The  read-write  statements  in  the  listing  are  written 
for  magnetic  tape  input  and  output  as  follows: 

(1)  Input  (READ  INPUT  TAPE  2) 

The  input  to  the  program  is  contained  on  one  card  (731C.3) 
and  includes  the  following  seven  quantities  in  order: 

I nput  quantity  Format 


SLNESS  (L/q>  ) 

E10.3 

ALPHA  (  c<  ) 

E10.3 

ECC  (e/q>) 

E10»3 

BETA  (  qb  ) 

£10.3 

A  (n) 

E10.3 

SIGY  (3y) 

E1C.3 

DELT  (AT) 

E10.3 

If  Subroutine  GRAPH  is  used  axis  scaling  and  multiple 
curve  data,  graph  title  data  and  curve  label  data  in 
addition  to  the  above  are  also  needed  in  the  input. 

(2)  Output  (WRITE  OUTPUT  TAPE  4) 

The  output  Includes  the  input  quantities  and  the  vari¬ 
able  quantities  for  every  fifteenth  point  as  follow: 
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Output  Quantity 

Format 

3LNESS  (lA>  ) 

F5o0 

ALPHA  (  ) 

F5.2 

ECC  ( e/ ) 

F5.3 

BETA  (  ) 

F4  *2 

A  (n) 

F5d 

SIGY  (Sy) 

F6  03 

DELT  (AT) 

F5..3 

T 

F5  02 

DISPL  (x/xt) 

El404 

STRN1  (€*) 

El404 

STRN2  (£.) 

E1404 

ST  RSI  (S,  ) 

E1404 

3TRS2  (Sz) 

214  A 

DEFL  (y/^  ) 

E14C4 

FORCE  (P/Pc  ) 

El  4  04 
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PROGRAM  COL  1 


C  PROGRAM  TC  DETERMINE  INELASTIC  COLUMN  RESPONSE  TO  HALF-SINE 
C  RAM  DISPLACEMENT  PULSE  APPLIEC  AT  END  WITH  SMALL  ECCENTRICITY 
C  TO  INITIALLY  STRAIGHT  COLUMN.  WHEN  RAM  CAUSES  TENSION  COLUMN 
C  EXHIBITS  FREE  VIBRATION. 

OD I  MENS  ION  T  (  900  ) , C I SPL ( 90 C) , C EFL ( 900 ) , STRN 1 ( 900 ) ,STRN2 ( 9C0), 
1  STRS1  ( 90C  )  ,STRS2(  900  ,  X<  4  )  ,  QCV  (  4 )  ,  EPS  II  (  4  )  ,  EP  S 12  (4  *,SIG1  (4)  , 
2 S I G2 ( 4  )  ,AK( 4,2) ,C (L )f FORCE (900) 

COMMON  A  ,  S I G Y 

READ  INPUT  TAPE  2 , 5  ,  SLNE SS , ALPH A , ECC , BET  A, A , S IGY , D ELT 
S  F0RMAT(7E1C.3) 

T  (  1 ) =C . C 
OEFL (  I )=0 .0 
VEL=0 . C 
STRN  1  (0=0.0 
STRN2  (0=0.0 
STRS1  (C  )  =0.0 
STRS2  (0=0. C 
C(  1  )  =  C.O 
C( 2)=C. 5 
C  (  3 )  =  C . 5 
C  (  4 )  =  1  .  C 
KTENS= 1 
KON  T  ACT= 1 
M  =  0 


l'/-5 


1  =  1 

40  CO  70  J  =  1 f  4 

TC  =  T (  I 1+C( J  )«DELT 

VC  =  DE  FL (  I )+C( J)*AK( J-l, 1 ) 

IF( 0. 5*BET  A  «TC-3.  14 1 59  3) SO, 55, 55 
50  X( J)=ALPHA*SINF(0.5*BETA*TC) 

GO  TO  60 
55  X(J)=C.O 

600EPSI 1  C J)  =  12. 1761/ (  SLNESS*«2  ) * ( X ( J ) -0. 5* C ECC** 2 )* VC **2* 
1  EtC*VC*(  1 .0-ECC)) 

I F ( J-  1)62,62,63 

62  CALL  STRESS(EPStl (J), STRN 1(1-1), STRS1 (1-1), SIGI(J) ) 

STRN1  ( I ) =E  PS  I  1 (1) 

STRS1 (  I  )  =S I G 1  (1 ) 

GO  TO  64 

63  CALL  STRESS(EPSI 1  ( J ) ,STRN1 (  l ) ,STRS1 ( I ) , SIG  1  ( J  )  ) 
640EPSI2(J)=12. 1761/ (  SLNESS**2  )*(X(J)-0.5*(ECC*«2)*VC**2- 

1ECC*VC*( 1 • 0+ECC ) ) 

IF(J-  1)65,65,66 

65  CALL  STRESS(EPSI2(J),STRN2( I-1) »STRS2( I— 1 ) , S IG2 ( J) ) 

STRN2 (  I  )  =  EPSI 2(1) 

STRS2 ( I  )  =  S I G2 ( 1 ) 

GO  TO  661 

66  CALL  STRESS( EPSI2 ( J) ,STRN2( I ) ,STRS2( I ),SIG2( J ) ) 

661  IFtKTENS-  1  )67 ,665 ,665 
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665  IFIK0NTACT-1 167  »6  94  »  694 


67  I  F ( M-  1  1675,68,685 

i 

675  XI =X( J ) 

P1  =  SIC-1  ( J  )  +S IG2 ( J  ) 

X( J)=CISPL{ 1-1 ) 

M=M+1 
GO  TC  60 

68  X2*X(J) 

P2*S  I  G 1 ( J 1 +S I G2 1 J ) 

X( JJ*X2-i P2/IP2-P1) ) * { X2-X 1 ) 

P  =  M  +  1 
GO  TO  6C 

685  P=SIG1 C  J  1+SIG21 J) 

IF( P)686,688,688 

686  I F ( 10.0**{-8.0)-AeSF(P) )6e7, 693,693 

687  X1=X{J) 

P1=P 

X(J)=X2-(P2/(P  2-P 1 1  )■ (X2-X1 1 
GO  TO  60 

688  I F  C 10.0**1-8.01 -A BSF(P) 1689,693,693 

689  X2  =  X ( J ) 

P2=P 

X(  J )  =  X2- ( P2/C P2-P 1 1 )*(X2-X1 ) 

GO  TO  6C 
693  N-0 
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KTENS=1 


GO  TO  698 

694  IF ( SI G 1 ( J ) +SIG2 ( J )) 695 , 698, 698 

695  IF( C. 5«BETA*TC-3. 141593)6955,696,696 
6955  KTENS=0 

GO  TO  675 

696  KONTACT  =  C 
GO  TO  675 

6 98000  V ( J  )=C. 04 10  639# (  SLNESS**2  )  » (  (  S  IG1  (  J )  *SIG2  ( J )  )*  (  VC*  1 . 0  )  - 
1 (SIG1 ( J  )- S  IG2  ( J  ) ) /ECC) 

AK ( J , 1)=DELT*( VEL+CI J) •DELT»CCV( J) > 

70  AK( J,2  )=CELT*DDVI J) 

DISPLt I )=X( 1 ) 

FORCE!  I  )  =  .0410639«(  SLNESS»*2  )•( STRS1( I  )*STRS2 (  I)  ) 

VEL=VEL* ( AK{ 1 , 2 )♦ 2. 0 • AK { 2 , 2 ) *2. 0* AK ( 3 , 2 ) ♦ AK ( 4 , 2 )  )/6.0 
OEFL (  1  +  1 )=DEFL l I)  ♦ ( AK ( 1 , 1  )  *2. 0* AK ( 2 , 1  ) *2 .0*AK ( 3 ,  1 ) *AK ( 4 , 1  >)/6.0 
IF (  1-900)75,80,80 
75  TCI  +  1  )  =  T(  I  )  *DELT 
1=1*1 
GO  TO  4 C 

80  WRITE  OUTPUT  TAPE  4,82 
82  FORMAT ( 1 H 1 ) 

WRITE  OUTPUT  TAPE  4 , 85 , SLNE S S , ALPHA , E CC , B ET A , A , S IG Y , C ELT 

» 

850F0RMAT ( 14H  OUTPUT  OATA  //3X 1 2HSLENDERNESS-F 5 . 0 ,2 X6HA LPHA=F 5 . 2 , 

12X4HECC=F5.3,2X5H8ETA=F4.2,2X2HA=F5. 1 ,2X5HSIGY«F6. 3,2X5H0ELT =F5. 3 
//  ) 


2 
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WRITE  OUTPUT  TAPE  4,90 


9CCFORMA T ( 4 X 4HT I  ME, 4 X12HD IS PLACEMENT, 4X10H0E FLECTION,  4X8HSTRAIN  1  , 
16X8HSTRESS  1 , 6X8HSTRA IN  2.6X8HSTRESS  2, 6X5HFORCE// ) 

00  92  1=1 ,900 

92  DEFL  (  I )  =  ECC»DEFL(  I  ) 

DO  95  1*1, 900,15 

950WRITE  OUTPUT  TAPE  4 , 98, T ( I ) ,0 ISPL (I) , DEFL ( 1 1 , STRNl (  I ) , STRS  1  (  I )  , 
1  STRN2( I ) ,STRS2(  I) , FORCE!  I ) 

98  FORMAT (3XF6.2»7E1  4.4) 

CALL  GRAPH1900.T,  STRM.8  ) 

CALL  GRAPH(900, T, STRN2,8) 

CALL  GRAPH(9C0,T, STRS1,8  ) 

CALL  GRAPH(900,T, STRS2.8) 

CALL  GRAPH! 900, T, 0EFL,8) 

CALL  GRAPH(900,T, FORCE, 8) 

STOP 

END 

SUBROLTINE  ST  RE SS ( E PS  I , E PS  I  A , S I GA , S I G ) 

COMMON  A , S I  GY 
I F { EPSI  120,20,21 

20  SIG=SIGA4(EPSI-EPSIA) 

GO  TO  33 

21  IFl EPSI-EPSIAJ22, 22,23 

22  S1G=SIGA4(EPSI-EPSIA) 
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GO  TO  33 


23  SIGB=$IGA  +  (EPSI-*EPSIA) 

8=3.0/(7.0*(SIGY**(A-1.0) )) 

S IG 1 =EPS I 

« 

2U  SIG2  =  SIG1-(SIGUB*(SIG1**A)-EPSI ) / ( A*B* < S IG 1  *  * ( A-l .0) ) ♦ 1 .0 ) 
IF(ABSF(SIG2*(10.  0**{-6.0) ) >-ABSF ( S IG2-S IG 1 ) >29,30,30 

29  SIG1 =S IG2 
GO  TO  24 

30  IF (SIG2-SIG8)31,3 1,32 

31  S I G=S IG2 
GO  TO  33 

32  SIG=SIGB 

33  RETURN 
ENO 
END 
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d.  Notation  for  Program  COL  2 

In  the  notation  for  Program  COL  2  the  subscript  J  re¬ 
fers  to  quantities  within  the  Runge-Kutta  prediction  loop 
which  are  calculated  for  intermediate  times  within  the  time 
increment  AT.  The  variable  quantities  ending  in  B  corre¬ 
spond  to  a  time  T  -  AT  and  the  variable  quantities  ending 
in  A  correspond  to  a  time  T.  The  notation  for  Program  COL  2 
is  as  follows: 


Computer  program  notation  for  Program  COL  2 


Computer 

Program 

Term 

Corresponding 
terms  in  table 
of  symbols 

Computer 

Program 

Term 

Corresponding 
terms  in  table 
of  symbols 

A 

n 

SIG1( J) 

s, 

ALPHA 

cx 

SIG2( J) 

Sa 

BETA 

$ 

SLNESS 

L/^ 

DDV ( J) 

i « 

V 

STRN1A  or  B 

DEPLA  or  B 

V 

STRN2A  or  3 

DEFLP 

yr/<^ 

ST  RSI  A  or  B 

s, 

DELT 

AT 

STRS2A  or  3 

sz 

DISPLA  or  B 

x/x£ 

STRS1M 

Max  S,  (tensile 
or  compressive) 

ECC 

STRS2M 

Max  Sz  (tensile 
or  compressive 

EPSI1( J) 

6, 

T 

T 

EPSI2( J) 

VELA  or  B 

V 

PM  AX 

X(J) 

x/xE 

SI  GY 

Sy 
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e.  Operating  instructions  for  Program  COL  2 

To  operate  Program  COL  2  use  the  cards  Including  Sub¬ 
routine  STRESS  in  the  order  shown  in  the  listing  followed  by 
one  blank  card,  then  the  input  data  card. 

The  read-write  statements  in  the  listing  are  written  for 
magnetic  tape  Input  and  output  as  follows: 

(1)  Input  (READ  INPUT  TAPE  2) 

The  input  to  the  program  is  contained  on  one  card  (7E10.3) 
and  includes  the  following  seven  quantities  in  order: 

Input  quantity  Format 

SLNESS(LA>  )  E10.3 

ALPHA  (  )  E10.3 

ECC  (e/^)  E10.3 


BETA  (  Q  ) 
A  (n) 


E10.3 


E10.3 


SIGY  (Sy) 
DELT  (AT) 


£10.3 


E10.3 


(2)  Output  (WRITE  OUTPUT  TAPE  4) 


The  output  includes  the  input  quatities  and  one  value 
of  each  of  the  variable  quantities  including  the  time 


T  of  occurrence  of  the  value  as  follows: 


Output  quantity 

SLNE3S  (L/^) 


Format 


F5.0 


ALPHA  ( cx.  ) 
ECC  (e/^  ) 
BETA  (  §  ) 

A  (n) 


F5.2 


E14.4 


F4.2 


F5.1 
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Output  quantity 

Format 

SIGY  (Sy) 

F6.3 

DELT  (AT) 

F5.3 

PMAX  (Pm/Pt) 

T1 

E14.4 

F6.2 

STRS1M  (SlfnaJ 

T2 

E14.4 

F6.2 

STRS2M  (SlmJ 

T3 

E14.4 

F6.2 

DEFLP  (yr/o) 

T4  ' 

E14.4 

F6.2 
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PROGRAM  COL  2 


C  PROGRAM  TO  DETERMINE  RESIDUAL  DEFLECTION  AND  MAXIMUM  FORCE 
C  IN  INELASTIC  COLUMN  RESULTING  FROM  HALF-SINE  RAM  DISPLACEMENT 
C  PULSE  APPLIED  AT  END  WITH  SMALL  ECCENTRICITY  TC  INITIALLY 
C  STRAIGHT  COLUMN.  WHEN  RAM  CAUSES  TENSION  COLUMN  EXHIBITS  FREE 
C  VIBRATION. 

ODIMENSICN  XU)  „OD  V  ( 4  ) ,  EPS  II C  4  )  ,  EPS  121  4  )  ,  S  IG 1  (  4 )  ,  SI  G2  (  4  > ,  AK  ( 4  ,2  )  , 
1CU> 

COMMON  A , S I GY 

READ  INPUT  TAPE  2  , 5  ,  SLNESS,  AL  PHA  ,  ECC,  BET  A ,  A  ,  S  I C-Y  ,D  ELT 
5  FORMATI7EI0.3) 

1C  T=0.0 

Cl SPL  A=0 . 0 
VEL  A=0 .0 
DEFLA=C.C 
STRN1 e-o.o 
STRN2  E=0 . 0 
STRS1 e=c.c 
STRS2B=0.0 
PMAX=C . 0 
STRS1 M=C . C 
STRS2  M  =  C  •  0 

i 

STRS2E=0 . 0 
WORK  =  C  .0 
C(1 )=C.C 
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C ( 2 ) =C  .5 


C(3)*0.5 
C  (  4  )  =  1.0 
KENERG®  1 
KT  ENS=  1 
KONT  ACT  =  1 
M= 0 

40  DO  70  J=1  ,4 
TC=T+C ( J ) *DELT 
VC=DEFLA*Cl J)*AK( J-1,  1  ) 

IF(0.5*BETA*TC-3. 141593)50,55,55 

50  X( J)=ALPhA*SINF(0.5*8ETA*TC) 

GO  TO  60 

55  X(J)=C.C 

600EPSI 1 (J)=12.1761/ (  SLNESS**2  > * ( X( J )-0. 5* ( ECC**2 )• VC* *2* 
1ECC*VC*( 1 .O-ECC)) 

I F (J— 1 )62,62,63 

62  CALL  STRESS(EPSI1  ( J ) , STRN IB, STRS 1 B,SI  G1 ( J ) ) 

STRN1 A  =  EPSI 1  (  1 ) 

STRS1 A=SIG1< 1 ) 

GO  TO  64 

63  CALL  STRESStEPSIl  (  J  ) ,  STRN  1  A,  STRS  1  A,  SI  GU  J)  ) 

640EPSI2 l  J )  = 12. 1761/ (  SLNESS**2  )* ( X ( J )-0. 5* ( ECC**2 )* VC* *2- 
.  1  ECC  *VC* ( 1 .O  +  ECC)) 

IF  (  J-  1  ).65 , 65, 66 
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65  CALL  STRESS(EPSI2 ( J) , STRN2B, STRS2B,SIG2< J)  ) 


STRN2  A  =  E  PS  I  2 (  1  ) 

STRS2  A  =  S I G2 l 1 ) 

GO  TO  661 

66  CALL  STRESS<EPSI2(J),STRN2A,STRS2A,SIG2(J>> 
661  IFIKTENS-l 167,665 ,665 

665  IF(KOKTACT— 1 167  ,694,694 

67  IF(M-1)675,68,685 
675  X1=X(J) 

P1=SIG1(J)+SIG2CJ  ) 

X( J)=CISPLB 
M  =  M+1 
GO  TO  6C 

68  X2=X(J) 

P 2  =  S I  C- 1  (  J1+SIG21J  ) 

XIJ) =X2- ( P2/( P2-P 1 ) ) «(X2-X1 ) 

M=M+1 
GO  TO  60 

685  P  =  SIG1  (  J)  +SIG2U) 

I F ( P ) 6  86  *  688 ,  688 

686  I F (  1 0 .  C*»  1-8.0) -A 8SFi P) } 6  87,693,695 

687  Xl=X(J) 

P1  =  P 

XU)  =  X2-(P2/1P2-P  i  )  1  *  (  X2-X 1 J 
GO  TO  60 
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688  I F  (  10.0«*(-*8.0)-A8SF(P)  )689,693»693 

689  X2=X(J) 

P2  =  P 

X(  J )  =  X2- ( P2/( P2-P 1 ) ) * ( X2-X1 ) 

GO  TO  6C 

693  M=0 
KTENS= 1 
GO  TO  698 

694  IF(SIG1(J)+SIG2(J ) >695,698,698 

695  IF (0. 5*BETA*TC-3. 141593)6955,696,696 
6955  KTENS=0 

GO  TO  675 

696  KONT ACT  =  C 
GO  TO  675 

6980DDV<J)=0. 04 10639*  (  SLNESS**2  )•( (SIG1  (J)+S«IG2tJ>  1MVC+1.C) 

I 

1 (SIG1 (J)-SIG2( J)> /ECC> 

AK( J,  1 )=CELT* ( VEL  A  +  C ( J)»DELT*0DV( J  )  ) 

70  AK(J,2)=DELT*DDVl J) 

IF (0.5*BETA*T-3.1 41593)702,704,704 
7020^0RK=fc0RK+0. 5*( < S TRS 1 A+S TRS2 A ) ♦ ( STRS 1 B+S TRS2B  )  )* 

1 (ALPHA*(SINF(0.5*BETA*T)-SINF(0.5*BETA*(T-OELT)  ) )) 

704  I F ( PMAX-C. 0410639* (  SLNESS**2  ) *  I STRS 1 A  +  STR S2 A ) ) 7 1 , 7 2 , 7 2 

71  PMAX^C .04  10639* (  SLNESS**2  ) * { STRS 1 A+STRS2A) 

T  1  =  T 

72  I F ( SIRS) V-STRS 1A)  73,74,74 
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73 

74 

75 

751 

752 

755 

756 

757 

758 

759 

76 

77 

78 

79 

80 


STRSlf*  =  STRS1A 


STRS1 E=STRS1M 
T2  =  T 

IF(ABSF(STRS2M)-ABSF{STRS2A) >75,751,751 
STRS2N  =  ST  RS2A 
T3  =  T 

I F( STRS2E-STRS2A) 752,755,755 
STRS2  E  =  ST  RS2A 
IF(KENERG-1 >756,76,76 
IF(VEIB)758»  759,7  57 
IF(VELA)759,7  59,8.1 
IF(VELA)81,759,759 
B=3.0/(7.0*( S IGY* • ( A- 1 .0 ) l ) 

ENERG 1=  ( A*B/ (A  +  l .0) ) *STR$ 1 E** ( A+ 1 . 0) *0 . 5*S7R  S  1 B  **  2 
ENERG2=  ( A*8 / ( A  +1  .0) ) *ST RS2E* * ( A  + 1 .0 ) +0. 5*STRS2B**2 
ENERGT=0. 0821278* (  SLNESS**2  > *< ENERG 1 +ENERG2 ) 

GO  TO  84 

IF ( KONTACT- 1 >77 ,8 1 ,81 
BENDM1=STRS1B-STRS28 
BEN0M2  =  STRS1 A-STR  S2 A 
I  F ( STRS1B  >78,78,79 
IF(STRS1A)81,80,80 
IF  (  STRS 1A  >80, 80,8  1 

DEFLP= ( DEFL A  -BEN CN2 / ( 8ENDM2- BENDM 1 ) * ( OEFL A-OEF L B ) )*ECC 
T4  =  T 
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KENERG=0 


81  T=T+DELT 

0  I SPLE=D I SPLA 
DISPLA=X(1) 

VELB=VEL A 

VEL A= VELB+ ( AK ( 1 ,2)+2.0«AK(2,2)+2.0*AK(3,2)+AK(4,2)  >/6.0 
DEFLB=0EFLA 

DEFLA  =  DEFLB  ♦( AMI , 1 ) +2 .0*AKC 2, 1 )*2. 0*AK ( 3, 1 ) +AK { 4, 1 ) ) / 6  .  C 

STRN 1 B  =  STRN 1 A 

STRN2  EsSTRN2A 

STRS 1 E=STRS1 A 

STRS2B=STRS2A 

GO  TO  40 

84  WRITE  OUTPUT  TAPE  4,85,SLNESS,  tCC  , BETA ,A , $1  GY ,0E LT 
850F0RMAT ( / // 1 4H  OUTPUT  OATA  //12HSLEN0ERNESS«F5.0,2X6H  ECC=F5.2, 
12X5HBETA=F4.2,2X2HA=F5. 1 ,2X5H$IGY*F6. 3, 2X5HOELT*F5 . 3/ / ) 

WRITE  OUTPUT  TAPE  4,86 

860 F0RMATI6X5H ALPHA, 7X1 1  HOE FORMAT  ION, 3X9HMAX  FORCF, 5X 1 2HMAX  STRESS  1, 
12X12HNAX  STRESS  2//) 

87  WRITE  OUTPUT  TAPE  4 ,88, ALPHA, OEFLP , PM AX, STRS1 Mr$TRS2M 

88  F0RMATI5E 14.4) 

WRITE  OUTPUT  TAPE  4 , 885 , T4, T1 , T2, T3  . 

885  FORMAT (20X1HIF6. 2, 1 H ) , 6X1 H ( F6 .2 , 1H ) ,6X1 H ( F6 . 2 , 1H ) , 6X 1 H ( F6. 2, 1H ) ) 

* 

WRITE  OUTPUT  TAPE  4 , 887, WCRK, ENERGT 
887  FORMAT { 25 X5HW0RK* E 1 2 .4, 1 0X7H ENERGY* E 1 2.4/) 
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STOP 


ENO 

SUBROUTINE  STRESS (EPS  I, EPSIA.S IGA, SIG) 

COMMON  A, SIGY 
IF(EPSI)20»20»21 

20  SIG=S IGA+(EPSI-EPSIA> 

GO  TO  33 

21  IF(EPSI-EPSIA)22, 22,23 

22  SIG=SIGA+(EPSI-EPSIA) 

GO  TO  33 

23  SIGB=SIGA+(EPSI-EPSIA) 

B=3.0/(7.0*(SIGY*»(A-1 .0) ) ) 

S IG1 =EPS I 

24  SIG2  =  SIGI-(SIGl+B*(SIG1**A)-EPS,n/(A«0«(Siqi*«(A-1  .01  +  1.0) 
IF(ABSF(SIG2*( 10. 0«* (-8.0)) )-ABSF( SIG2-SIG1 1)29,30,30 

29  SI G1=SIG2 
GO  TO  24 

30  IF(SIG2-SIG8)31,31,32 

31  S I G=S IG2 
GO  TO  33 

32  S IG*S IG8 

33  RETURN  .  . 

END 

END 
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APPENDIX  V 


Maximum  Force.  Residual  Deflection  and  Maximum  Stress  Data 

Tabulated  on  the  following  pages  are  the  maximum  force, 
residual  deflection,  and  maximum  stress  data  for  the  three 
slenderness  ratios  and  two  eccentricity  ratios  for  all  values 
of  CK  and  Q  which  were  Investigated  for  pulse  loading.  The 
values  of  maximum  lateral  deflection  that  were  available  are 
also  Included  in  the  tabulation.  These  values  were  only  avail¬ 
able  from  the  solutions  to  computer  Program  COL  1  used  to  ob¬ 
tain  graphical  results.  Thus,  only  a  limited  number  of  these 
values  are  Included  in  the  tabulation. 

The  dimensionless  time  Increments  AT  used  for  the  numer¬ 
ical  integration  were  the  same  for  a  given  value  of  in 
all  the  solutions  for  both  computer  programs.  These  values 
are  as  follows: 


$ 

AT 

0.01 

0.10 

0.05 

0.10 

0.10 

0.08 

0.20 

0.04 

0.50 

0.02 

1.00 

0.01 

The  values  of  AT  are  more  than  adequate  to  maintain  the 
accuracy  of  the  numerical  scheme.  The  minimum  number  of 
increments  for  a  natural  period  of  first  mode  lateral  vibra¬ 
tion  of  the  column  in  the  numerical  process  is  62.  The  values  of 
AT  were  selected  to  display  the  complete  history  of  the 
dynamic  response  to  the  pulse  loading  for  the  graphs  shown 
in  App.  Ill  using  the  maximum  900  points  permitted  in  Sub- 
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routine  GRAPH  * 

The  values  of  maximum  stress  tabulated  on  the  following 
pages  are  maximum  compressive  stresses  if  there  Is  no  sign 
before  the  number,  or  maximum  tensile  stresses  if  there  Is 
a  negative  sigh  before  the  number.  In  either  case  they  are 
the  maximum  absolute  values  of  stress  that  occurred  In  Spring 
1  or  Spring  2  of  the  column  model. 
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Maximum  Force,  Residual  Deflection  and  Maximum  Stress  Data 


2  I/^-  1Q0  e/^«  0,05 


$ 

Maximum 

Force 

(p„/p6  ) 

Residual 

Deflection 

(y*-/?) 

Maximum 

Deflection 

(yn/^) 

Stress  1 
Maximum 

Stress  2 
Maximum 
< 

0.01 

0.953 

2.12x10"* 

0.00282 

0.00074 

0.05 

1.001 

2.24x10”* 

—  • 

0.00283 

0.00078 

0.10 

1.121 

2.50x10'* 

— 

0.00287 

0.00100 

0.20 

1.319 

2.77x10-* 

— 

0.00289 

-0.00129 

0.50 

1.65$ 

2.40x10"* 

0.00286 

-0.00233 

1.00 

1.935 

1.69x10"* 

— -  ' 

0.0027$ 

0.00213 

<=<-  2.5 


e/^a  0.05 


$ 

Maximum 
Force 
<*U/fV ) 

Residual 

Deflection 

Maximum 

Deflection 

(y»n/£) 

Stress  1 
Maximum 

Stress  2 
■Maximum 
(s2.max) 

0.01 

0.953 

5.89x10"* 

~ 

0.00312 

0.00084 

0.05 

1.029 

6.45x10"* 

—  - 

0.00315 

-0.00099 

0.10 

1.189 

7.66x10"* 

— 

0.00320 

-0.00132 

0.20 

1.470 

9.41x10"* 

— 

0.00327 

-0.00188 

0.50 

1.958 

8.88x10”* 

_  — i 

0.00325 

ro. 00275 

1.00 

2.331 

9.00x10"*’ 

0.00329 

-0.00263 
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Maximum  Force,  Residual  Deflection  and  Maximum  Stress  Data 


«<-  3  .  '  L/^>-  100  e/^»  0.05 


$ 

Maximum 
Force 
(Pm/Pc  ) 

Residual 
Deflection 
(yr/^)  . 

Maximum 

Deflection 

(ym/9) 

Stress  1 
Maximum 

(^irnax) 

Stress  2 
Maximum 

0.01 

0.954 

1.14x10"’ 

--- 

0.00334 

-0.00110 

0.05 

1.069 

1.28x10"' 

-  —  — 

0.00337 

-0.00130 

0.10 

1.252 

1.54x10-' 

—  —  — 

0.00344 

-0.00167 

0.20 

1.602 

1.98x10“ 

—  — _ 

0.00352 

-0.00236 

0.50 

2.237- 

2.20x10“ 

— —  » 

0.00356 

-0.00310 

1.00 

2.633 

2.05x10“ 

—  : 

0.00360 

-0.00313 

«*>-  3.5 


e/^a  0.05 


> 

Maximum 

Force 

(p„/r, ) 

Residual 

Deflection 

(yr/^) 

Maximum 

Deflection 

(y»n/^) 

Stress  1 
Maximum 

Stress  2 
-Maximum 

( ax) 

0.01 

0.955 

1.81x10“ 

0.00349 

-0.00131 

0.05 

1.110 

2.05x10“ 

0.00354 

-0.00155 

0.10 

1.314 

2.47x10“ 

_  — _ 

0.00360 

-0.00195 

0.20 

1.724 

3.25x10“ 

—  —  - 

0.00370 

-0.00268 

0.50 

2.447 

3.83x10“ 

0.00377 

-0.00340 

1.00 

2.851 

2.57x10“ 

0.00380 

-0.00352 
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Maximum  Force,  Residual  Deflection  and  Maximum  Stress  Data 


c<-  4  I/^-100  e/^„  0.05 


$ 

Maximum 

Force 

(pVpe  ) 

Residual 

Deflection 

(yr/9) 

Maximum 

Deflection 

(ym/9) 

Stress  1 
Maximum 

(SUnnax) 

Stress  2 
Maximum 

0.01 

0.956 

2.54x10"' 

■p— ..  „ 

0.00361 

-0.00149 

0.05 

1.147 

2.39x10“' 

— 

0.00366 

-0.00175 

0.10 

1.362 

3.49x10“' 

0.00373 

-0.00219 

0.20 

1.336 

4.56x10"' 

— 

0.00333 

-0.00231 

0.50 

2.617 

5.57x10“' 

— 

0.00392 

-0.00367 

1.00 

3.011 

2.93x10“' 

0.00394 

-0.00331 

o<»  1  70  e/^a  0.05 


Maximum 

Force 

(p„/p, ) 

Residual 

Deflection 

(yr/^) 

Maximum 

Deflection 

(y»n/^) 

Stress  1 
Maximum 
( 

Stress  2 
-Maximum 
( 6  2,  max) 

0.01 

0.341 

3. 17x10" 1 

- • 

0.00315 

0.00149 

0.05 

0.340 

3 . 19xl0-i 

— 

0.00315 

0.00150 

0.10 

0.347 

3. 13x10“ z 

— 

0.00315 

0.00153 

0.20 

0.363 

4. 13x10“ 1 

— 

0.00324 

0.00163 

0.50 

0.934 

.2.33x10"* 

0.00312 

0.00190 

1.00 

0.932 

6.24x10"’ 

— 

0.00272 

0.00224 
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Maximum  Force,  Residual  Deflection  and  Maximum  Stress  Data 


1.5 


70 


e/^.  0.05 


$ 

Maximum 
Force 
(Pm/P6  ) 

Residual 

Deflection 

(y*-/Q)  • 

Maximum 

Deflection 

Stress  1 
Maximum 

Stress  2 
Maximum 

( szma^ 

0.01 

0.841 

2.10x10“' 

- - 

0.00381 

0.00149 

0.05 

0.854 

2.09x10"' 

— 

0.00381 

0.00153 

0.10 

0.898 

2.10x10“' 

0.00381 

0.00164 

0.20 

0.987 

2.10x10“' 

- - 

0.00381 

0.00184 

0.50 

1.206 

2.04x10“' 

— 

0.00381 

-0.00260 

1.00 

1.337 

1.03x10“' 

- —  ■ 

0.00364 

0.00314 

2  Iy^  70  e/^=  0.05 


Maximum 

Force 

) 

Residual 
Deflection 
(y  r/^>) 

Maximum 

Deflection 

(y*n/0 

Stress  1 
Maximum 

Stress  2 
•Maximum 
(s2.max) 

0.01 

0.841 

3.79x10“' 

- — 

0.00404 

0.00149 

0.05 

0.878 

3.86x10“' 

*  * 

0.00405 

0.00.158 

0.10 

0.954 

3.96x10“' 

— 

0.00406 

0.00176 

0.20 

1.080 

4.05x10“' 

- . 

0.00407 

-0.00228 

0.50 

1.377 

3.43x10“' 

0.00403 

-0.00368 

1.00 

1.516 

2.22x10“' 

— 

0.00402 

0.00364 
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Maximum  Force,  Residual  Deflection  and  Maximum  Stress  Data 


c<-  2.5 


L/^-  70 


e/^„  0.05 


§ 

Maximum 
Force  1 

(p~,/pt ) 

Residual  . 
Deflection 
(y*-/?) 

Maximum 

Deflection 

(ym/^) 

Stress  1 
Maximum 

Stress  2 
Maximum 

0.01 

0.641 

5.36x10“' 

— 

0.00416 

0.00149 

0.05 

0.697 

3.52x10“'  . 

— 

0.00419 

0.00161 

0.10 

0.995 

5.77x10“' 

- - 

0.00421 

-0.00202 

0.20 

1.167 

.  5.99x10“' 

1.67 

0.00423 

-o„ 00306 

0.50 

1.469 

5.03x10*' 

— 

0.00422 

-0.00413 

1.00 

1.619 

3.06x10“' 

- • 

0.00422 

0.00393 

1  ,  50  e/^a  0.05 


Maximum 

Force 

UU/r. ) 

Residual 

Deflection 

(yr/^). 

Maximum 

Deflection 

(y»n/0 

Stress  1 
Maximum 
( 

Stress  2 
■Maximum 

( SlLmax) 

0.01 

0.565 

2.67x10“' 

— 

0.00417 

0.00234 

0.05 

0.592 

2.69x10“' 

— 

0.00416 

0.00240 

0.10 

0.614 

2.70x10h 

— 

0.00419 

0.00251 

0.20 

0.655 

2.63x10“' 

— 

0.00417 

0.00275 

0.50 

0.731  • 

2.52x10“' 

— 

0.00419 

0.00329 

1.00 

0.774 

6.60x10“* 

0.00396 

0.00366 
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Maximum  Force,  Residual  Deflection  and  Maximum  Stress  Data 
c<-  1.5  1/^m  50  e/^»  0.05 


$ 

Maximum 

Force 

(p,yps ) 

Residual 

Deflection 

(y»-/Q) 

Maximum 

Deflection 

(ym/^) 

Stress  1 
Maximum 

( ^imax) 

Stress  2 
Maximum 

0.01 

0.535 

5.17x10*' 

— 

0.00446 

0.00234 

0.05 

0.612 

5 . 16x10“' 

0.00446 

0.00251 

0.10 

0.654 

5. 13x10“' 

0.00446 

0.00275 

0.20 

0.708 

5.05x10”' 

- -  . 

0.00446 

0.00307 

0.50 

0.822 

4.03x10”' 

— 

0.00445 

0.00383 

1.00 

0.864 

1.60x10“' 

— - 

0.00435 

0.00413 

o<»  2  1^5-  50  e/^»  0.05 


Maximum 
Force 
(?„/!>«  ) 

Residual 

Deflection 

(yr/^) 

Maximum 

Deflection 

(y*n/£) 

Stress  1 
Maximum 

( ®imax) 

Stress  2 
Maximum 
( ^  2.  max) 

0.01 

0.535 

7.27x10"' 

— 

0.00461 

0.00235 

0.05 

0.626 

7.34x10"' 

— 

0.00462 

0.00257 

0.10 

0.679 

7.37x10”' 

— 

0.00462 

0.00287 

0.20 

0.756 

7.18x10”' 

— 

0.00462 

0.00336 

0.50 

0.873 

4.93x10”' 

•  0.00461 

-0.00454 

1.00 

0.913 

2.07x10”' 

— 

0.00456 

0.00439 
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Maximum  Force,  Residual  Deflection  and  Maximum  Stress  Data 


c<-  2.5  I/^-  50  e/^«  o.05 


$ 

Maximum 
Force 
(Pm/Pc  ) 

Residual 

Deflection 

Maximum 

Deflection 

(y.n/9) 

Stress  1 
Maximum 

(^irnax) 

Stress  2 
Maximum 

0.01 

0.586 

9.14x10”' 

— 

0.00472 

0.00236 

0.05 

0.645 

9.30x10”' 

— 

0.00473 

0.00269 

0.10 

0.700 

9.40x10“' 

— 

0.00473 

0.00300 

0.20 

0.794 

8.98x10“' 

— 

0.00474 

-0.00413 

0.50 

0.909 

6.17x10”' 

— 

’  0.00475 

-0.00475 

1.00 

0.948 

2.41x10"' 

0.00471 

0.00456 

c<-  2  Iy^-  100  e/^a  0.01 


Maximum 

Force 

(p„/p« ) 

Residual 

Deflection 

(yr/^) 

Maximum 

Deflection 

(y»n/£) 

Stress  1 
Maximum 
( ) 

Stress  2 
•Maximum 

( ^  2.  max) 

0.01 

1.003 

2. 5 6x10” 1 

— 

0.00237 

0-.  00100 

0.05 

1.192 

3.03x10"* 

— 

0.00292 

0.00113 

0.10 

1.376 

3. 2 9x10“ 1 

1.64 

0.00295 

0.00149 

0.20 

.  1.616 

3.30x10"* 

1.81 

0.00295 

-0.00180 

0.50 

1.911 

3.07x10"*- 

1.93 

0.00294 

-0.00231 

1.00 

1.988  . 

•  2. 95x10" 3 

0.69 

0.00250 

0.00235 
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Maximum  Force,  Residual  Deflection  and  Maximum  Stress  Data 


2.5  •  I/^-  100  0.01 


$ 

Maximum 
Force 
(Pm/Ps  ) 

Residual 

Deflection 

(yr/9). 

Maximum 

Deflection 

Stress  1 
Maximum 

Stress  2 
Maximum ' 

0.01 

1.017 

6.93x10“* 

— 

0.00317 

0.00100 

0.05 

1.258 

8.89x10"* 

— 

0.00325 

-0.00139 

0.10 

1.486 

1.04x10"' 

— 

0.00331 

-0.00185 

0.20 

1.813 

1.17x10"' 

--- 

0.00334 

-0.00250 

0.50 

2.264 

1.11x10"' 

—  —  — 

0.00335 

-0.00286 

1.00 

2.423 

2.29x10"* 

0.00306 

0.00287 

c<«  3 


e/^a  0.01 


Maximum 

Force 

(p„/p5 ) 

Residual 

Deflection 

(Yr/^) 

Maximum 

Deflection 

(y*n/0 

Stress  1 
Maximum 
( 

Stress  2 
■  Maximum 
( 2  2.  max) 

0.01 

I.O32 

1.30x10"' 

- - 

0.00338 

-0.00115 

0.05 

1.320 

1.72x10"' 

- -  ’ 

0.00348 

-0.00173 

0.10 

1.578 

2.07x10"' 

- - 

0.00354 

-0.00224 

0.20 

1.971 

2.47x10"' 

— 

0.00360 

-0.00274 

0.50 

2.535 

2.11x10h 

0.00359 

-0.00326 

1.00 

2.744 

7.27x10"*- 

-— 

0.00346 

0.00327 

V-10  .  • 


Maximum  Force, 

* 

<=*-  3.5 


$ 

Maximum 
Force  ' 

(pVps  ) 

Residual 

Deflection 

(yr/9) 

Maximum 

Deflection 

(ym/9) 

Stress  1 
Maximum 

( Sur*ax) 

Stress  2 
Maximum 

( ^irnax^ 

0.01 

1.04S 

2. 00x10"  ‘ 

— 

0.00353 

-0.00135 

0.05 

i.3ao 

2.69x10"' 

— 

0.00363 

-0.00199 

0.10 

1.666 

3.26x10“' 

: - 

0.00370 

-0.00255 

0.20 

2.106 

3.a6xlO“' 

— 

0.00377 

-0.00292 

0.50 

2.739 

3.00x10“' 

- — 

0.00376 

-0.00356 

1.00 

2.961 

1.35x10“' 

—  . 

0. 00372 

0.00354 

o-C.  •  4  100  e/^a  0  •  01 


?■ 

Maximum 

Force 

( V*. ) 

Residual 

Deflection 

(yv/^) 

Maximum 

Deflection 

Stress  1 
Maximum 

Stress  2 
-  Maximum 

( S^max) 

0.01 

1.062 

2.74x10“' 

--- 

0.00364 

-0.00151 

0.05 

1.433 

3.70x10“' 

— 

0.00375 

-0.00220 

0.10 

1.751 

4.52x10“' 

2.67 

0.00333 

-0.00273 

0.20 

2.230 

5.49x10“' 

2.97 

0.00390 

-0.00303 

0.50 

2. 394 

3.94x10“' 

3.54 

0.0033a 

-0.00333 

1.00 

3.114 

1.93x10“' 

2.70 

0.00390 

0.00374 

APPENDIX  V 

Residual  Deflection  and  Maximum  Stress  Data 

l/<£-  100  e/^»  0.01 
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Maximum  Force,  Residual  Deflection  and  Maximum  Stress  Data 


1  70  0,01 


$ 

Maximum 
Force 
<P~/Pt  ) . 

Residual 

Deflection 

(yr/9) 

Maximum 

Deflection 

(y.n/9) 

Stress  1 
Maximum 

Stress  2 
Maximum 

0.01 

0.916 

2.54xl0‘* 

— 

0.00309 

0.00197 

0.05 

0.927 

2. 89x10- z 

— 

0.00313 

0.00201 

0.10 

0.947 

3.63x10”* 

0.49 

0.00320 

0.00206 

0.20 

0.973 

2.28x10”* 

0.45 

0.00306 

0.00217 

0.50 

0.991 

3. 94x1 0”3 

0.20 

0.00264 

0.00232 

1.00 

0.994 

Iil3xl0"3 

0.08 

0.00252 

0.00242 

<=<-  1.5 


L/^  70 


a/^»  0.01 


Maximum 

Force 

(ruA* ) 

Residual 

Deflection 

(yr/^) 

Maximum 

Deflection 

(y*n/^) 

Stress  1 
Maximum 

(^irnax) 

Stress  2 
. Maximum 
(Sz.max) 

0.01 

0.922 

2.24x10”' 

— 

0.00383 

0.00199 

0.05 

0.994 

2.24x10”' 

— 

0.00383 

0.00212 

0.10 

1.072 

2.22x10"' 

— 

0.00383 

0.00230 

0.20 

1.181 

2.17x10“' 

— 

0.00383 

0.00261 

0.50 

1.333 

1.60x10”' 

- - 

0.00377 

0.00315 

1.00 

1.364 

2.11x10”* 

— 

0.00345 

0.00334 
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Maximum  Force,  Residual  Deflection  and  Maximum  Stress  Data 


2 


1/^-  70 


e/^Ja  0.01 


$ 

Maximum 

Force 

(Pm/Ps  ) 

Residual 

Deflection 

(yr/^) 

Maximum 

Deflection 

(ym/9) 

Stress  1 
Maximum 

( sirnax) 

Stress  2 
Maximum 

0.01 

0.930 

4.01x10“' 

mm  mm  mm 

0.00406 

0.00200 

0.05 

1.041 

4.17x10"' 

— 

0.00408 

0.00222 

0.10 

1.151 

4.17x10“' 

— 

0.00408 

0.00250 

0.20 

1.296 

4.02x10“' 

— 

0.00408 

-0.00309 

0.50 

1.497 

3.04x10“' 

—  • 

0.00408 

0.00359 

1.00 

1.542 

5.20xl0“i 

— 

0.00388 

0.00380 

c<-  2.5  Iv^«  7Q  e/^»  0.01 


Maximum 

Force 

(Pm/re ) 

Residual 

Deflection 

(y  v-/^) 

Maximum 

Deflection 

(y*«/£) 

Stress  1 
Maximum 

Stress  2 
-Maximum 

( ^  z.  max) 

0.01 

0.939 

5.62x10“' 

0.00420 

0.00201 

0.05 

1.079 

6.02x10“' 

• 

0.00423 

0.00231 

0.10 

1.210 

6.17x10“' 

1.85 

0.00424 

-0.00292 

0.20 

1.362 

6.03x10“' 

1.95 

0.00425 

-0.00364 

0.50 

1.594 

3.95x10“' 

1.99 

0.00425 

-0.00396 

1.00 

1.642 

7. 93x10  “2- 

0.72 

0.00413 

0.00405 
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Maximum  Force,  Residual  Deflection  and  Maximum  Stress  Data 


o<.-  1 


L/^  ■  50 


e/^a  0.01 


§ 

Maxi nun 
Force 

(?m/?S  ) 

Residual 

Deflection 

(yr/9) 

Maximum 

Deflection 

(*»/?> 

Stress  1 
Maxi nun 

( ^lrnax^ 

Stress  2 
Maxi nun 

I8*™*) 

0.01 

0.633 

2.71x10"' 

— 

0.00419 

0.00234 

0.05 

0.658 

2.74x10"' 

— 

0.00420 

0.00297 

0.10 

0.689 

2.75x10“' 

0.67 

0.00421 

0.00313 

0.20 

0.730 

2.65x10“' 

0.75 

0.00421 

0.00337 

0.50 

0.775 

1.15x10“' 

0.48 

.0.00403 

0.00367 

1.00 

0.782 

1.78xlO“l 

0.10 

0.00384 

0.00378 

o  1,5  Iy^»  50  e/^s  0,01 


Maxi nun 
Force 

(?„/?, ) 

Residual 

Deflection 

( >v/^>) 

Maximus 

Deflection 

(y»n/ 

Stress  1 
Maximum 

( ^irna.x^ 

Stress  2 
. Maximum 
( S^max) 

0.01 

0.636 

5.30x10"' 

— 

0.00447 

0.00286 

0.05 

0.637 

5.22x10“' 

— 

0.00447 

0.00311 

0.10 

0.735 

5.02x10“' 

1.12 

0.00447 

0.00338 

0.20 

0.791 

4.57x10“' 

1.10 

0.00447 

0.00369 

0.50 

0.859 

2.62x10"' 

0.90 

•  0.00443 

0.00411 

1.00 

0.870 

3 . 66x10“* 

0.16 

0.00426 

0.00422 
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Maximum  Force,  Residual  Deflection  and  Maximum  Stress  Data 

■  1  L/^  m  50  6/^  a  0*01 


$ 

Maximum 

Force 

ou/p« ) 

Residual 

Deflection 

(y*~/<?) 

' Maximum 
Deflection 

(Ynr,/^) 

Stress  1 
Maximum 

Stress  2 
Maximum 

0.01 

0.633 

2.71x10“' 

— — 

0.00419 

0.00284 

0.05 

0.658 

2.74x10“' 

0.00420 

0.00297 

0.10 

0.689 

2.75x10“' 

0.67 

0.00421 

0.00313 

0.20 

0.730 

•  2.65x10“' 

0.75  • 

0.00421 

0.00337 

0.50 

0.775 

1.15x10-' 

0.48 

.0.00403 

0.00367 

1.00 

0.782 

1*  78x10“ * 

0.10 

0.00384 

0.00378 

c<-  1.5 


e/^a  0.01 


Maximum 

Force 

uu/p* ) 

Residual 

Deflection 

(yr/^) 

Maximum 

Deflection 

(y*n/£) 

Stress  1 
Maximum 

Stress  2 
-Maximum 
( s2.max) 

0.01 

0.636 

5.30x10“' 

mm  mm 

0.00447 

0.00286 

0.05 

0.687 

5. 2 2x10“ 1 

mm  mm  mm 

0.00447 

0.00311 

0.10 

0.735 

5.02x10-' 

1.12 

0.00447 

0.00338 

0.20 

0.791 

4.57x10“' 

1.10 

0.00447 

0.00369 

0.50 

0.859 

2.62x10"' 

0.90 

•  0.00443 

0.00411 

1.00 

0.870 

3 • 66x10“* 

0.16 

0.00426 

0.00422 
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Maximum  Force,  Residual  Deflection  and  Maximum  Stress  Data 


OC-  2  L/^>-  50  °*01 


$ 

Maximum 

Force 

(Pm/Ps  ) 

Residual j 
Deflection 
(yr/^) 

Maximum 

Deflection 

(y.n/<$) 

Stress  1 
Maximum 
( S  \  max  ) 

Stress  2 
Maximum 

0.01 

0.641 

7.44x10"' 

— 

0.00463 

0.00269 

0.05 

0.709 

7.44x10"' 

— 

0.00463 

0.00323 

0.10 

0.766 

7.17x10"' 

— 

0.00464 

0.00354 

0.20 

0.632 

6.27x10“' 

— 

0.00463 

0.00391 

0.50 

0.907 

3.44x10"' 

0.00463 

0.00435 

1.00 

0.919 

4. 44x10" \ 

— 

0.00449 

0.00446 

c*.-  2.5  50  ®/^  0.01 


Maximum 

Force 

(p„/p« ) 

Residual 

Deflection 

(yr/^) 

Maximum 

Deflection 

(y*n/0 

Stress  1 

Maximum 

( 

Stress  2 
. Maximum 
( stmax) 

0.01 

0.646 

9.35x10"' 

— 

0.00473 

0.00291 

0.05 

0.729 

9.45x10"' 

- — 

0.00475 

0.00334 

0.10 

0.769  . 

9.19x10"' 

1.65 

0.00476 

0.00367 

0.20 

0.663 

7.94x10"' 

1.60 

0.00476 

-0.00427 

0.50 

0.941 

4.03x10"' 

1.26 

0.00476 

0.00452 

1.00 

0.952 

5.25xlO“L 

0.22 

0.00465 

0.00462 
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Evaluation  of  the  Approximations  and  of  the  Accuracy  of 
the  Numerical  Scheme 

From  a  consideration  of  the  maximum  lateral  deflections 
encountered  in  the  solutions  the  approximations  made  in  de¬ 
veloping  the  equations  for  the  column  model  may  be  evaluated. 
The  maximum  value  of©  for  all  solutions  was  less  than  0.08 
radian.  Based  on  a  maximum  ©  =  0o08  radian  the  errors  are 
as  follows: 

(a)  0.01  per  cent  error  in  (L/2-e  tan©  )  as  a  result 
of  neglecting  e  tan© 

(b)  0.2  per  cent  error  in  Eraine=e(l  -cos-©-  3  by  con¬ 
sidering  sin©=©  and  neglecting  e(  1-cos-©) 

(c)  0.25  per  cent  error  in  tan© and  0.11  per  cent 
error  in  sin©  as  a  result  of  considering 
tan  ©  =  sin©=© 

(d)  less  than  0.06  per  cent  error  in  (1-cos©)  as  a 
result  of  considering  l-cos©=£©4 

From  this  the  approximations  which  have  been  made  in  deriv¬ 
ing  the  equations  for  the  column  model  are  reasonable  and 
produce  insignificant  error  in  the  results. 

As  an  additional  check  for  cumulative  error  a  compar¬ 
ison  of  the  ram  work  done  on  the  column  model  with  the  net 
energy  input  to  each  of  the  springs  (a  combination  of  strain 
energy  of  plastic  deformation  and  the  existing  elastic  strain 
energy)  is  made  in  all  computer  solutions  for  the  results 
shown  in  Figs.  10  thru  14.  These  agree  within  1  per  cent 
thus  verifying  the  compatibility  of  the  equations  and  c'neck- 
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ing  the  accuracy  of  the  numerical  scheme  for  the  time  Incre¬ 
ments  chosen.  Work-energy  considerations  are  not  Included 
In  the  computer  program  for  obtaining  graphical  results  but 
the  dimensionless  time  increment  AT  used  for  the  numerical 
integration  in  both  programs  for  a  given  value  of  was  the 
same. 
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